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X Abstract 
(D 

The electromagnetic form factors of the proton and the neutron in the timehke region are re- 
^ viewed. In view of the forthcoming experimental projects devoted to investigate these observables, 

^ we present the current status of the field and we emphasize the relevant role, that accurate mea- 

surements, refined phenomenological analyses, as well as microscopic models will play with the 
goal of achieving deeper insights into the structure of the nucleon and its inner dynamics. 
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1 Introduction 

J> The investigation of the nucleon structure through electromagnetic (em) probes has played a central 
*^ role in our present understanding of strong interactions. In particular, the measurements of em nucleon 
form factors (FF's) in the spacelike (SL) region, which started with the pioneering work of Hofstadter 
[1], have given the indisputable evidence of the non-elementary nature of the proton. The possibility, 
offered by SL photons, of achieving an accurate description of the three-dimensional structure of the 
nucleon has motivated many and successful experimental efforts, which resulted in an accurate data 
base of the SL nucleon FF's (see Ref. P] for a list of recent reviews). With the advent of a new 
experimental technique - the so-called polarization-transfer method - a striking fall-off of the proton 
electric FF relative to the magnetic one, which largely follows a dipole behavior, has been discovered. 
This relatively recent achievement impressively demonstrates the relevance of the polarization degrees 
of freedom in the studies of the nucleon structure. 

The transition to the timelike (TL) region, widening in such a way the investigation of the nucleon 
FF's, opens unique possibilities for studying peculiar features of the nucleon state. In particular, in the 
TL region, one can exploit the hadronic components of the virtual photon state, that become physical, 
according to the proper energy thresholds. A very popular model of those hadronic components is the 
vector meson dominance model (see Ref. [3] for a recent review), that allows i) to develop a simple 
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picture of the production mechanism of a physical hadron-antihadron pair, and ii) to accurately identify 
quantities, which are containing relevant details on the nucleon/antinucleon states. The information 
accessed in the TL region could appear less intuitive than the electric and magnetic distribution densities 
of the nucleon, but indeed they play a relevant role for shedding light on the long-range behavior of 
the strong interactions. TL FF's are indeed more sensitive to the discrete hadron spectrum (given the 
pole structure that occurs) and to the presence of many transition amplitudes of physical processes. 
For instance, though the SL nucleon FF's can be described by models, that take into account only 
few vector mesons, the pole pattern shown by the experimental TL FF's cannot be reproduced without 
considering all the relevant vector mesons. To better understand the appealing feature of the TL region, 
let us go back to the fifties, when the theoretical studies on this issue began. In particular, in a paper 
by Federbush, Goldberger and Treiman [4] , it was pointed out that a bridge can be established between 
the studies of the nucleon em FF's and the transition amplitudes of physical processes, like em decays of 
vector mesons (VM's) as well as the production of hadronic states with well-defined quantum numbers. 
This link was found through the dispersion relations that allow to express the FF's in terms of an 
integral over their imaginary parts, that are non-vanishing only in the TL region. At that time, the 
mesonic and nucleonic degrees of freedom were the relevant ones in the field of the strong interactions, 
and therefore the mentioned em VM decays and hadronic transition amplitudes were considered as the 
relevant inputs for a phenomenological analysis. Our modern paradigm of strong interactions. Quantum 
ChromoDynamics (QCD), allows us to understand those physical quantities, through a microscopical 
description in terms of quark and gluon degrees of freedom. It should be noticed that microscopic 
models could represent an effective tool for evaluating the above mentioned quantities (see, e.g., the 
analyses in Refs. O El E]), since Lattice QCD is greatly challenged in the TL region, differently from 
the SL case where electric and magnetic density distributions begin to be addressed [8j. 

Summarizing what has been illustrated above, it is clear that a measurement of the imaginary part 
of the TL nucleon FF's, that can be achieved by exploiting the polarization degrees of freedom, is 
highly desirable, with a final goal of constructing a more detailed database of the nucleon FF's, that 
does contain not only information on the electric and magnetic distribution densities, but also on other 
physical hadronic processes. 

At present, the experimental knowledge of the TL nucleon FF's unfortunately does not yet allow 
to achieve this goal. Compared to the SL sector, where precision measurements of FF's have been 
achieved on the percent level, the data base of TL FF's is rather scarce. Recent measurements of the 
BaBar collaboration [9] have somewhat improved the knowledge of the proton FF. The BaBar result 
concerning the ratio of the electric to the magnetic FF is however still largely limited by statistical 
uncertainties, and it is furthermore in conflict with a previous measurement at LEAR [TU]. This only 
demonstrates the need for a new generation of high-precision measurements of the em FF of the proton 
and neutron. Indeed, new facilities at Novosibirsk, Beijing, and FAIR/Darmstadt recently came into 
operation or are about to start in the coming years. From these facilities we can expect significantly 
improved results. 

The aim of the work is therefore to give an updated (maybe not complete) view to the present 
status of both the experimental and theoretical investigations in the field. Moreover, we hope that this 
review could represent an opportunity, allowing to share the common expertise, and to compare the 
different experimental as well as theoretical techniques, which have been worked out in the past. But, 
after all, the main goal is to show the wealth of information which is provided by the TL nucleon FF's 
for understanding in depth the non-perturbative regime of QCD. In view of this, we will emphasize 
that only measurements of the relative phases will fully determine the nucleon FF's in the TL region, 
so that the maximal phenomenological knowledge can be reached. Therefore, this review appears to be 
timely, given the near- and mid-term future experimental programs. 

The review is organized as follows. Sect. 2 is devoted to the experimental studies, in particular in 
Subsect. 2.1 an overview of the facilities as well as the experimental techniques is given. In Subsect. 2.2, 
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the present experimental data base is presented, and in Subsect. 2.3 the future facihties and perspectives 
are discussed. In Sect. 3 we present the theoretical background. We start in Subsect. 3.1 with a brief 
review of the SL em FF's; in Subsect. 3^ the general formalism for investigating the TL nucleon FF's 
is presented. It follows in Subsect. 3^ a discussion of the threshold energy region; in Subsect. 3^ the 
path how to obtain the dispersion relations between the real and the imaginary part of the nucleon FF's 
is shown; in Subsects. 3.5 - 3.9, several theoretical approaches are finally reviewed. In Sect. |4] some 
conclusions are drawn and future perspectives are presented. 
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2 The experimental investigation of TL nucleon form factors 



In the year 1972 the first measurement of a TL nucleon FF was performed at the e^e~ colhder ADONE 
in Frascati using the process e^e~ — )■ pp [llj. This historically first result was obtained with an optical 
spark chamber setup at a center-of-mass energy (cm. energy) of y/s = 2.1 GeV/c. In the following 
years a series of measurement campaigns were performed at the electron-positron colliders ADONE 
with the FENICE experiment [121 [13], as well as at the Orsay colliding beam facility (DCI) with the 
detectors DM1 [H] and DM2 [151 US]. The em FF of the proton was explored by these facilities from 
nearly production threshold {^/s = 2Mp ■ c, where Mp is the mass of the proton) up to cm. energies 
of 2.4 GeV/c. Precision measurements were also obtained with the BES-II experiment [17| at BEPC, 
and with CLEO [I8l at CESR. The information concerning the TL FF of the neutron is especially 
scarce, with only the FENICE experiment so far having been able to perform this measurement, cf. 
Refs. [ni[IS|. 

First attempts to measure the proton FF using the inverse reaction pp — )■ e^e~ date back to the mid 
1960's. Indeed, first upper limits stem from antiproton beam experiments at BNL [20] and CERN [21] . 
The discovery of this reaction was finally possible using an antiproton beam at PS/CERN in 1976 [22]- 
Antiproton experiments were later continued with great success at LEAR/CERN with the PS170 ex- 
periment [ini 1231 El] and at FNAL [23 ESI EZ] • 

In the e~^e~ experiments mentioned above, em FF's have been measured using the so-called energy 
scan method, i.e. by systematically varying the cm. energy of the e~^e~ collider. Around the beginning 
of the 21*^* century, e~^e~ particle factories came into operation, such as the B-factory PEP-II at SLAC, 
which was operated at a cm. energy corresponding to the mass of the T(4S) resonance of 10.6 GeV. 
It was realized that the use of events with photon radiation from the initial state (ISR) appears to 
be a copious source of hadronic final states with invariant masses below the actual cm. energy of the 
collider. Competitive results of the proton [0] as well as of the AA, AS", and E'^S^ FF's [2B] have been 
achieved at the BaBar experiment in the course of the years. The BaBar data set does not only feature 
the best statistical and systematic precision achieved to date, but it is also spanning the entire energy 
range of interest from threshold up to 4.5 GeV. The BELLE collaboration at the Japanese B-factory 
project KEK-B has been using the ISR-technique in the charmonium energy region and has measured 
the process e+e" A+A;: [23- 



2.1 Experimental techniques and facilities 



We report in this Subsect. on the various experimental techniques (see also Fig. (|T])) and facilities, as 



outlined above. A compilation of the experimental results follows in Subsect. 2.2 An outlook concerning 



future perspectives in the field will finally be discussed in Subsect. 2.3 



2.1.1 Energy scan experiments in e+e annihilation 



The standard technique for measurements of hadron production in electron-positron annihilation is the 
so-called energy scan, in which the cm. energy of the collider, s/s, is varied systematically, see Fig. ([ija). 
At each energy point a measurement of the associated cross section is carried out. In the case of nucleon 
pair production, e+e~ — )■ NN, the associated FF's are related to the total cross section in one-photon 



approximation (see Subsect. 3.2 for a more detailed discussion) according to: 



6S 



+ 



2Ml 



where a is the electromagnetic fine structure constant, /3 the nucleon velocity, Cn{s) is the s-wave 
Sommerfeld-Gamow factor, that takes into account the Coulomb effects at threshold (see Eq. (22) 
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and below for more details), and the nucleon mass. and are the electric and magnetic 
Sachs FF's, respectively, which depend on the square four-momentum transfer of the virtual photon, 
q^. Neglecting radiative corrections, the four-momentum transfer and the cm. energy of the collider 
are identical, cf. Eq. ([T]): 



= s. 



In none of the energy scan experiments performed so far, the Sachs FF's could be disentangled by means 
of a measurement of the differential cross section (see below), and the data are therefore shown in terms 
of an effective FF defined as follows: 



|Gfjtf)l 



3g2 



CNiq' 



2M% 



g2 + 2M 



where the rightmost-hand side shows the relation with the Sachs FF's in one-photon approximation. 

We want to stress that any experimental measurement of is independent of any assumption on 
\Ge\ or Q ^iid it quantitatively indicates how much the experimental cross section differs from a 



point-like one (cf. Eq. (27)). Summarizing, the cross section a^^+g^-^^fj can be written as a product of 
a point-like cross section and |G^(g^)p, that contains information on the transition amplitude from a 
virtual photon state to a NN state (i.e, in a naive picture, I7*) — )■ \qqqqqq) — )■ \NN), through all the 



allowed paths, see also Subsect. 3.2) 



The cross section itself is obtained by measuring the total number of nucleon pairs, N'nn, after 
background subtraction (Abkg) and by normalizing to the total integrated luminosity (£int)- Since 
detectors cannot cover the full solid angle, the geometrical acceptance (egeom) and detection efficiencies 
(cdetect), as Well as radiative corrections ((5rad) need to be considered: 

_^fNN-^fhkg 111 

(^e+e-^NN — ' ' ' , r • l-jj 

■'-int Cgeom ^detect -L ~r Orad 

Table ([T]) shows a compilation of nucleon FF measurements, which have been obtained via the energy 
scan technique. While the pp channel has been measured by several experiments (Refs. pT| [T2| [T3| [HI 
[T5| US]), data on the nn final state so far exist only from FENICE |T2l [13]. The experiment DM2 also 
measured at one single energy point the process AA [15]. As can be seen from Table ([!]), the energy 
range from production threshold up to ~ 3 GeV has been covered via energy scan experiments. Typical 
energy steps in the individual scanning campaigns vary from 40 MeV - in the case of the Frascati and 
Orsay experiments - up to 100 MeV in the case of BEPC |l7j. At 3.67 GeV, one single energy point has 
been published by CLEO at Cornell [18]. With integrated luminosities well below 1 pb~^ per scan point 
(in almost all cases) and given the small cross section for the process e+e~ — )■ pp, the collected statistics 
typically was very low. Even close to threshold, where the cross section for baryon-pair production 
is highest (but still below 1 nb), the statistics per scan point did not exceed 25 events. We briefiy 
discuss in the following the methodology used at the experiments DM1 / DM2 (Orsay), the FENICE 
experiment (Frascati), and the BES experiment (Beijing). 

• At the Orsay colliding beam facility DCI the nucleon FF measurements were performed with 
the detectors DM1 and DM2. A sketch of the DM2 detector is shown in Fig. (|2]). It already 
resembles the typical structure of modern collider experiments with large geometrical acceptance. 
DM2 consisted of several layers of cylindrical multiwire proportional chambers (MWPCs) inside 
a magnetic field of solenoidal shape. In two scan campaigns (four and six scanning points each) 
the energy range between 1.925 GeV and 2.25 GeV could be covered; in total, approximately 170 
pp events were selected. 



^ In literature, it is often stated that the experimental value of \G^Q{q^)\ is obtained assuming \Ge\ — \Gm\i which is 
incorrect. 



5 



Figure 1: Diagrammatic representation of the experimental processes used for the measurement of 
timehke nucleon FF's: (a) e+e^ annihilation scan experiments; (b) pp annihilation; (c) the Initial State 
Radiation technique at e~^e~ colliders; in all cases the form factor is measured as a function of the 
square four-momentum transfer of the virtual photon coupling to the baryon pair. 
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Figure 2: Detector DM2: (1) Proportional wire Figure 3: Detector FENICE in a view transversal 
chamber, (2) Drift chamber, (3) Cerenkov coun- to the beam direction. Figure taken from Ref. [13] . 
ters, (4) Scintillation counters, (5) Photon detec- 
tor, (6) Muon detector, (7) End cap detector. Fig- 
ure taken from Ref. [30]. 

• The non-magnetic FENICE detector (cf. Fig. ([3])) at the ADONE collider was designed with 
the primary goal to measure the neutron FF. The main components of the experiments were 
Limited-Streamer- Tube (LST) modules as tracking devices, scintillation counters as triggering and 
timing devices and thin iron plates as distributed converters for the antineutron detection [13] . 
The specific signature of an antineutron annihilation 'star' together with the long TOF between 
production and annihilation, were mainly used to detect anti-neutrons. No dedicated neutron 
detection was performed due to the relatively low neutron detection efficiency. The neutron FF 
has been measured from values close to threshold up to 2.44 GeV. A precise determination of the 
proton FF in the same energy range has been possible as well. 

• From an energy scan (10 scanning points) of the BEPC collider in the range 2.0 GeV/c < ^/s < 
3.07 GeV/c, the BES-II experiment could extract the proton FF, extending significantly the energy 
range covered before at Frascati and Orsay. BES-II was a typical multiple-purpose detector with 
large acceptance, consisting of a vertex chamber, a large cylindrical drift chamber, a TOF system 
as well as a photon detector setup. The whole apparatus was embedded in a solenoidal coil. The 
precision of the proton FF measurement at BEPC was mostly limited by statistics. In total, 80 
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Ji/Xp. 


rveaction 


Year 
Publ. 


Scan 
Points 


riange [ij-e v j 


-^int 
[pb-1] 


Ji/vents 


ivei. 


DM1 


e'^e~ —J- pp 


1979 


4 


1.925 -2.180 


0.4 


~ 70 


m 


DM2 


e~^e~ —J- pp 


1983 


6 


1.975 - 2.25 


0.5 


~ 100 


m 


DM2 


e~^e~ — )■ pp 


1990 


1 


2.4 


0.2 


7 


m 


DM2 


e+e" AA 


1990 


1 


2.4 


0.2 


4 


m 


ADONE 73 


e~^e~ — pp 


1973 


1 


2.1 


0.2 


25 


m 


FENICE 


e+e^ — 7- 


1993 


2 


2.0-2.1 


<0.1 


27 


m 


FENICE 


e^e^ — > 


1993 


1 


2.1 


0.1 


28 


m 


FENICE 


e+e^ —7- 


1994 


4 


1.9-2.4 


0.3 


70 


m 


FENICE 


e^e^ — )■ nfi 


1998 


5 


1.9 - 2.44 


0.4 


74 


m 


FENICE 


e+e" — )■ pp 


1998 


1 


2.1 


<0.1 


7 


m 


BES-II 


e+e^ —J- pp 


2005 


10 


2.0-3.07 


5 


80 


m 


CLEO 


e"'"e~ — )• pp 


2005 


1 


3.671 


21 


16 


m 



Table 1: Summary of experimental results obtained at e~^e colliders via the energy scan technique. 

e~^e~ — )■ pp events have been detected. The systematic uncertainty was close to 10% in all scan 
points. 



Shower det . _ LH2 tarQet 




Figure 4: Experiment PS 170 at the antiproton storage ring LEAR. Figure taken from Ref. [TO]. 



2.1.2 pp annihilation 

Starting in the 1960's, antiproton beams became available, which were also used to search for the 
process pp — )■ e^e~, see Fig. ([l|b). The antiproton beam is scattered on a hydrogen target in a fixed- 
target configuration. The momentum transfer squared, g^, of the em FF in this case is accessible by 
measuring the invariant mass of the lepton pair, or, alternatively, it is calculated for a given antiproton 
beam momentum assuming the target protons to be at rest. First upper limits for the TL proton FF 
have indeed been obtained at relatively high energies at BNL [20] and CERN pi] in ranges between 
5 (GeV/c)^ and 8 (GeV/c)^. For historical reasons we report FF results from e"'"e~ experiments as 
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Figure 5: Differential cross section for the reaction pp — )■ e~^e~. The angle 6* denotes the angle between 
the outgoing lepton and the antiproton beam direction in the cm. reference frame. The differential 
cross section has been measured by PS170 for five different antiproton momenta. The full line is a fit 
to the differential cross section assuming \G^^\ = the dashed line shows a fit to data with \G\^\ 

and |G|;/G^^| used as free parameters. Figure taken from Ref. pUj . 



a function of the cm. energy of the collider, ^/s, while results from pp annihilation experiments are 
quoted as a function of g^. 

In the year 1976 - 3 years after the first Frascati measurement of a TL nucleon FF in e~^e~ scattering 
- the first positive result was achieved by the Mulhouse-Strasbourg- Torino collaboration [22] with 
an antiproton beam at PS/CERN. In total, 25 candidate events, which were produced by stopped 
antiprotons in a liquid hydrogen target, as well as 5 additional events, in which the annihilation process 
took place 'in-fiight', had been recorded. The experimental setup consisted of several optical spark 
chambers surrounded by an array of scintillators. Due to the low antiproton beam momenta, for the 
first time the proton FF could be measured at values almost at the production threshold. 

Indeed, pp annihilation experiments allow for a complementary access to the TL proton FF. Dif- 
ferently to e~^e~ colliders, the production threshold does not feature a limited phase space for the 
outgoing particles, and furthermore, one does not have to deal with the experimental issues related to 
low-momenta protons in the target region. On the other hand, no baryon FF's different from protons 
can be accessed in pp annihilation experiments, and the rejection of background with hadronic final 
states is an experimental challenge. 

With the advent of a dedicated storage ring for low-momenta antiprotons at CERN, the LEAR 
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ring, significant progress could be achieved also for measurements of the TL proton FF. A dedicated 
experiment PS170 was able extract the FF at 9 different values from threshold up to 4.2 (GeV/c)^ 1101 
l23l Elj. Not only the beam intensity was increased with respect to the previous CERN experiment, 
but also the instrumentation was more advanced. The PS170 detector consisted of several layers of 
MWPCs as tracking detectors, Cerenkov detectors, and hodoscopes inside a 1.2 T magnetic field. It 
was possible to separate signal events from two-body processes involving hadrons. Fig. (|4]) shows a 
sketch of the detector. In a sophisticated analysis, the yield of the e^e~ final state was normalized to 
the yield of the two-body hadronic reactions pp — tt^tt^ and pp — )■ K^K~ . The cross section for the 
process pp — )■ e^e~ was then computed using the known ratio of hadronic two-body final states to the 
total pp annihilation cross section, cf. Ref. [10] for more details. 

For antiprotons stopped in the liquid hydrogen target - i.e. for events very close to threshold - a large 
statistics of almost 2, 000 events could be collected, exceeding the previously available statistics by more 
than one order of magnitude. The extracted FF was however limited by systematic uncertainties. Due 
to the fact that the stopped antiproton is forming a pp atom, only the configuration with J^^ = 1 
quantum numbers is relevant for the em FF. Hence, the fraction of L = 1 quantum states, as well 
as the ratio of the spin singlet to triplet probabilities need to be known. Overall, an uncertainty of 
approximately 10% was achieved at threshold. 

The large angular acceptance together with the relatively high statistics did also allow for a first 
measurement of the differential cross section using 'in- flight' events. In one-photon approximation, this 
differential cross section is given by 



da na^ ^ . 2 
: Cp[q 



dcosO 4/^1 



4M^ 

\GW)\' (1 + cos^^) + \Gl{q')\\m'9 



(4) 



where 6 is the angle between the proton and the beam direction in the cm. system, and k the antiproton 
momentum in the CM system ( see also Eq. ( [2T| . From the measurement of the differential cross section 
(cf. Fig. ([5])) and using Eq. Q, the ratio IG^/G^/I could be extracted. The experiment PS170 has 
measured this ratio for five different values, although with relatively large statistical and systematic 
uncertainties of ~ 30%. The result for the ratio will be presented in the next Subsect.. 

In the 1990's and at the beginning of the new century, measurements of the TL proton FF have 
been obtained also at FNAL, where a high-intensity, stochastically-cooled antiproton beam used to 
be in operation. Using a hydrogen gas target, the experiments E761 and E835 covered the high-g^ 
range between 8.2 (GeV/c)^ and 18.2 (GeV/c)^. Those were the highest momentum transfers achieved 
at the time, exceeding the measurements at BES-II in electron-positron annihilation. At the highest 
energy point, = 18.22 (GeV/c)^, only an upper limit could be established. In contrast to the low- 
momenta antiproton experiments performed at LEAR, at FNAL the cross section measurement could 
be normalized to the integrated luminosity, which was monitored by a dedicated monitoring system 
using elastically scattered antiprotons as a reference process. The cross section is determined using 
Eq. pi) - with the only difference that a different kinematic flux factor has to be considered - , and 
Eq. ([sf with A/^vtv to be replaced with AAee, where A/'ee corresponds to the number of pp — )■ e+e~ events. 

Table ^ shows a summary of all results from pp annihilation experiments for extracting TL nucleon 
FF's. 



2.1.3 Initial State Radiation 

Around the year 2000, a new generation of high-luminosity electron-positron colliders came into oper- 
ation, which were explicitly designed to operate at fixed cm. energies, ^/s, corresponding to the mass 
of either the 0(lO2O)-resonance - in the case of the 0-factory DA$NE in Frascati - or the mass of the 
T(4S')-resonance - in the case of the B-factories PEP-II at SLAG and KEK-B in Tsukuba. As those 
resonances provide the highest statistics for coherent production of pairs of K- and B-mesons, the main 
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JTVcdC LlUll 


Year 
Publ. 


Scan 
Points 




JZj VcllLo 




M.S.T. Coll. 


pp — )■ e~^e~ 


1976/77 


2 


near threshold 


34 


m 


PS170 


pp — )■ e~^e~ 


1991 


4 


near threshold 


~ 2000 


[23] 


PS170 


pp — 7- e+e~ 


1991 


4 


1.94-2.05 


~ 1300 


m 


PS 170 


— )■ e"'"e" 


1994 


9 


threshold — 2.05 


~ 2000 


m 


E760 


pp — )■ e"'"e~ 


1993 


3 


3.0-3.6 


29 


m 


E835 


pp — e+e~ 


1999 


4 


3.0-3.8 


144 


m 


E835 


pp — )■ e"^e~ 


2003 


2 


3.4-3.5 


66 


m 



Table 2: Summary of experimental results using the reaction pp — )■ e~^e . 

physics motivation for the above mentioned facilities was the measurement of CP-violation in K- and 
B-meson systems. 

Energy scans over wide energy ranges were technically not possible at the first generation 0- and 
B-factories. It was realized, however (see Refs. [311 |32l [33] ) , that the very high luminosities did allow for 
a complementary approach to the standard energy scan. Events with photon radiation from the initial 
state (ISR), cf. Fig. ([ijc), lead to a reduction of the invariant mass of the virtual photon, s', and hence 
allow for a measurement of the hadronic FF's at energies below ^/s. By measuring the invariant mass 
of the hadronic system, Mhadr, the entire energy range below the cm. energy of the collider becomes 
accessible: y/s' < y^. By neglecting effects of final state radiation (FSR), the following relation holds 
for the four- momentum transfer q"^: 

= ^hadr ■ = S'. 

Comprehensive reviews of the ISR method, which is also called Radiative Return technique, can be 
found in Refs. [31] and [35]. The method was both used for measurements of hadronic reactions with 
mesons in the final state, as for the measurement of baryon pair production: 

e+e" BB'jisK. 

The BaBar experiment at the SLAC B-factory PEP-II has measured the TL proton FF below 4.5 GeV, 
and has furthermore investigated the following three hyperon final states: AA, AS°, and E^S'^. Before 
BaBar no experiment had the sensitivity to measure the latter two channels. The full BaBar data set 
comprises an integrated luminosity of ~ 500 fb"^, collected between 1999 and 2008. BaBar results 
of mesonic final states, i.e. of final states containing two pions, three pions, four pions etc. are very 
valuable input for evaluating the hadronic contributions to {g — 2)^ and the running fine structure 
constant Aahad(^I)- Notably, a measurement of the TL pion FF from threshold up to 3 GeV with 
< 1% precision in the p peak region was achieved [36]. Furthermore, a series of interesting structures 
have been observed in the mass spectra, e.g. a destructive interference effect in the 6tt final states [37] 
at an energy being very close to the pp threshold. Among the major results of the ISR programme at 
BaBar, we want to mention also the discoveries of the y(4260) [SB] and y(2175) [SH] resonances. 

On the theoretical side, the ISR process was calculated within QED up to Next-to-Leading Order 
(NLO) and the probability for ISR photon emission can be expressed by means of a radiator function 
Hs{s') [321 ESI HQ]. Notice that high energies of the ISR photon correspond to small values of s'. 
Neglecting again effects of FSR (see Ref. [H] concerning this issue), i.e. assuming Mhadr ■ c = y/s', the 
non-radiative cross section o"hadr = cr(e^e~ — )■ hadrons) can be extracted from the measured radiative 
cross section crhadr+7 = <7{e~^e~ — > hadrons + 7isr) by using this radiator function: 

^^c/g'^^ ^ ^hadr(s') ■ Hs{s') ■ ^. (5) 
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Figure 6: ISR luminosity in the case of BaBar kinematics for L\ 
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The calculation of the radiator function is usually performed with the Monte Carlo method; a dedicated 
Monte Carlo generator for ISR processes has been developed within the PHOKHARA package, which 
simulates a large number of hadronic reactions, including also the baryonic final states pp7, nn^y^ 
and AA7 |13]. PHOKHARA simulates the full NLO-ISR radiative corrections [B]. Within the BaBar 
collaboration, besides PHOKHARA, also the AFKQED Monte Carlo package has been developed, which 
is simulating FSR corrections by means of PHOTOS [15] and which uses the structure function approach 
for higher order ISR corrections. The AFKQED generator is based on an early version of PHOKHARA, 
called EVA, cf. Refs. [321 SS]. 

Due to the specific energy dependence of the ISR photon, the radiator function Hs{s') decreases 
steeply with s'. In Fig. (|6]), the so-called ISR luminosity is shown, which is the product of the radiator 
function times the integrated luminosity: AsrI^') = Hs{s') ■ £int- The specific values for BaBar {y/s = 
10.6 GeV, £int = 500 fb~^) have been used in Fig. ([6|. The geometrical acceptance (egeom) and detection 
efficiency (edetect) have been considered as well (e = egeom ■ Cdetect)- 

The ISR luminosity allows to compare the number of ISR events produced for a given final state and 
for a given mass value s' with the direct production in an energy scan; £isr depends of course on the 
chosen bin size and decreases with finer binning. At the pp production threshold, the ISR luminosity 
is found to be 1.7 pb~^ (assuming a bin width of 100 MeV/c^) and therefore exceeds the statistics 
collected at ADONE and DCI. Apart from the large statistics, the ISR method offers some additional 
advantages over conventional energy scan measurements: 

• A very valuable feature is the fact that in one single experiment the entire energy range of interest 
can be covered. This avoids the notorious problems associated with combinations of different 
data sets, which are usually having different normalization uncertainties. At BaBar, the proton 
FF could be measured over a very wide mass range from production threshold up to y/s' = 4.5 
GeV/c, which has never been possible before and which exceeds the energy range covered by 
BES-II and the FNAL experiments. 

• In the baryon FF measurements performed at BaBar, the ISR photon has been tagged in the em 
calorimeter (EMC) at relatively large polar angles. Events with small angle ISR photons, which a 
priori have a much higher cross section, lead to an event signature, in which the hadronic system is 
emitted opposite to the ISR photon. The hadronic particles are hence flying outside the detector 
acceptance. The lower the invariant of the hadronic system is, the higher is the probability to lose 



11 




Figure 7: BaBar measurement of the process 
e^e~ — )■ pp-y. Dependence of the detection effi- 
ciency on cos 6p, where dp is the proton angle mea- 
sured in the pp rest frame with respect to the ISR 
photon direction. Right: Dependence of the Mpp 
mass resolution on Mpp. Figure reprinted with per- 
mission from Ref. [35], ©2011 American Physical 
Society. 



Figure 8: BaBar measurement of the process 
e+e~ —J- pp'j. Dependence of the Mpp mass res- 
olution on Mpp. Figure reprinted with permission 
from Ref. |35], ©2011 American Physical Society. 



at least one of the hadronic particles in the final state. However, for the tagged photon approach 
this back-to-back signature has its own merits, namely it leads to the fact that the geometrical 
acceptance for pp events is very fiat compared to usual energy scan experiments. This can be seen 
in Fig. ([t]), where the angular acceptance for e -|- e~ — )■ pp^ events is shown. The almost uniform 
angular distribution reduces the model dependence of the cross section measurement due to the 
unknown \G^^/G\.^\ ratio, which represents a major limitation in the data analysis of energy scan 
experiments. 

• Furthermore, for ISR events - differently from standard e~^e~ experiments - the momenta of the 
final state particles at threshold are non-zero. This does not only reduce the systematic effects, 
which are usually unavoidable when dealing with low-momenta particles, but it leads also to a 
flat reconstruction efficiency as a function of Mhadr = Mpp, which is especially important if the FF 
at threshold is investigated. In Fig. ^ the mass resolution as a function of the hadronic mass is 
shown. At low masses, the mass resolution is as small as 1 MeV/c^. In order to improve the mass 
resolution and to reject background channels, a kinematic flt is usually performed. 

We conclude that the ISR technique is a very competitive method for the measurement of baryon TL 
FF's. Moreover, the high statistics and accuracy of the BaBar measurement of the proton FF [H] allowed 
a measurement of the \G^^/G\.j\ ratio in flve energy bins. Like in the case of the PS170 experiment, 
the angular distributions are investigated to disentangle and However, when using the ISR 

technique, the existence of an additional high-energetic photon in the flnal state needs to be taken 
into account. In the BaBar analysis, the differential cross section with respect to the cosine of the 
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Figure 9: ISR measurement of the process e'^e — t- pp with the BaBar detector. The cos 6 distributions 
are presented for different pp mass regions: (a) 1.877-1.950 GeV/c^ (b) 1.950-2.025 GeV/c^ (c) 2.025- 
2.100 GeV/c2, (d) 2.100-2.200 GeV/c^, (e) 2.200-2.400 GeV/c^, (f) 2.400-3.000 GeV/c^. The points 
with error bars show data distributions after background subtraction. The histograms are fit results: the 
dashed histograms show the contributions corresponding to the magnetic form factor; the dash-dotted 
histograms show the contributions from the electric form factor. Figure reprinted with permission from 
Ref. [35], ©2006 American Physical Society. 
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Figure 10: The BaBar detector at the B-factory PEP-II at SLAG: (SVT) Sihcon Vertex Tracker, (DCH) 
Drift Chamber, (EMC) Electromagnetic Calorimeter, (DIRC) Detector of internally reflected Cerenkov 
Light, (IFR) Instrumented Flux Return. 



angle has been fitted (cf. Fig. ([9])), where 6p is the angle between the proton in the pp rest frame 
and the momentum of the pp system in the e~^e~ rest frame (remember that PEP-II is an asymmetric 
collider built for CP studies in the B- meson sector). In Ref. [13] it has been suggested to use a different 
approach, namely to consider the proton angular distribution in the hadronic rest frame, with the z-axis 
aligned with the direction of the ISR photon and the y-axis in the plane determined by the beam and 
the photon directions. 



Fig. (10) shows a sketch of the BaBar detector. Charged particles are detected in the BaBar 
tracking system, which comprises a five-layer silicon vertex tracker (SVT) and a 40-layers cylindrical 
drift chamber (DCH) operating in a 1.5 T solenoidal magnetic field. The momentum resolution for a 1 
GeV/c track is (Tpjpt = 0.5%. Particle identification is provided by a an internally reflected Cerenkov 
detector (DIRC), as well as by dE/dx information from the tracking system. The energy and position of 
photons is measured by a CsI(Tl) em calorimeter (EMC); muons are identified by a dedicated detection 
system inside the iron coil (instrumented flux return, IFR). No attempts have been made at BaBar to 
measure the detection efficiency of neutrons inside the EMC; only the proton-antiproton final state as 
well as hyperon channels, which decay into protons and pions, have been investigated. 



Table (|3j) gives a summary of the ISR measurements of TL baryon FF's, which have been performed 
so far. Besides the BaBar results reported above, the BELLE experiment at the Japanese B-factory 
project KEK-B has investigated the process e^e~ — )■ A+A~ in the energy range between 4.5 GeV and 
5.4 GeV using the ISR technique [29] . The integrated luminosity collected at BELLE is exceeding 1, 000 
fb~^ and is hence twice as large as the luminosity available at BaBar. Until now, the ISR programme at 
BELLE has been entirely devoted to the production of charm and charmonium resonances in the mass 
range above 4 GeV/c^. In the case of the A^A~ analysis, a resonance at M = (4634]t7(stat.)l8(syst.)) 
MeV/c^ and a width of F = (92t^°(stat.)t^?(syst.)) MeV has indeed been found. 
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Figure 11: Cross section measurements of the 
process e~^e~ — )■ pp obtained at e~^e~ colliders: 



Figure 12: The proton effective form factor 
extracted from measurements at e'^e'^ colliders 



BaBar [9], FENICE [13], DM2 [151 [16], DM1 [H], (Refs. as in Fig. ((TTj)) and at antiproton accel- 
ADONE73 [n], BES [H], CLEO [IS]. erators: PS170 [lOj, E835 i26ii27j, E760 [25]. 



2.2 Experimental results 

In this Subsect. the experimental results obtained in in the extraction of TL baryon FF's are sum- 
marized. While the main focus of this review is the physics of the TL FF of the proton and of the 
neutron, we also present the results obtained by BaBar on hyperon FF's. Furthermore, we discuss some 
phenomenological aspects related to the FF measurements. The in-depth discussion of the theoretical 
background follows in Sect. [3| 



2.2.1 TL proton electromagnetic form factor 

In Fig. ( pTj ) the various measurements of the cross section, performed at e^e~ colliders, are summarized, 
see Tables and g. The energy range between the production threshold and 4.5 GeV has been 
covered. Fig. ( |l2| ) shows the effective FF for those e"'"e~ data sets and includes in addition the results 
obtained via proton-antiproton annihilation at CERN and FNAL, cf. Table (|2]). We have omitted the 
older CERN results of Refs. [23l [Ml [22] which are superseded by the high-statistics measurements of 
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Table 3: Summary of TL baryon FF results obtained via the ISR technique. 
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Figure 13: Proton effective form factor measurements obtained in the energy range close to the pp 
threshold region: BaBar [9], PS170 [10], FENICE [I3], DM1 [H]. 



experiment PS 170 (Ref. [TO]). In general we can identify a good agreement among the different data 
sets. Nevertheless, the complex shape of the TL proton FF is largely not understood and has lead to 
many speculations. We briefly summarize the most relevant observations. 

• The FF shows a very steep rise toward threshold, which can be clearly identified in the two data 



sets of BaBar and PS170. A zoom to the threshold region is shown in Fig. 13 (linear vertical 
scale), in which the BaBar spectrum is plotted in a finer binning to take account of the rapid FF 
change. Although the BaBar and PS170 spectra seem to indicate some normalization issues, the 
threshold enhancement is clearly visible also in this plot. It is interesting to note that a similar 
behavior can be seen also in the pp invariant mass spectra of very different physics processes, 
e.g. B-decays (S° -> D^*^pp, 5+ K^pp) measured at BELLE [IT] and BaBar [48j as well 
as the decay J/ip ^ pp'j measured at BES |19]. It has been speculated whether the threshold 
enhancement might be due to the existence of a hypothetical, narrow resonance with a mass 
just below threshold. Indeed, such a narrow NN bound state could give rise to a dip in the 
energy dependence of the cross section e+e~ — )■ hadrons as a result of an interference effect with 
a broad resonance (SO]. Several experiments had observed such a dip in the energy spectrum of 
the Gn state [51j and more recently also BaBar has confirmed this observation in the 3(7r+7r^) and 
2(7r+7r^7r'^) exclusive states [37]. We will discuss the physics of the threshold region in more detail 
in Subsect. 13.31 



Fig. (14) shows the energy dependence of Ggff the energy range from threshold to 3.01 GeV 



(upper plot), and from 2.58 GeV to 4.5 GeV (lower plot). The events from J/ip and 'ip{2S) decays 
to pp are subtracted from the contents of the corresponding bins. Two rapid decreases of the FF 
and cross section near 2.25 GeV and 3 GeV are seen by BaBar. Rosner has pointed out [52j that 
these steps are just below the respective thresholds for pA(1232) and A^(1520)iV(1520) systems, 
respectively, and suggests an s-wave threshold effect to be responsible for these structures. 

Finally, a comparison of the asymptotic behavior of the TL proton FF at high q"^ with the corre- 
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Figure 14: Proton effective FF measurements plotted in linear scale for the energy range 1.85 — 3.01 
GeV (upper plot) and 2.58 - 4.5 GeV (lower plot): BaBar [9], FENICE [13], DM2 [ISl HE], DM1 [H], 
ADONE73 [IT], BES [H], CLEO [18], PS170 [ID], E835 [261 EI], E760 [25]. 



sponding SL FF, measured in e~p scattering, represents another puzzle. While perturbative QCD 
calculations and the application of the Phragmen-Lindeloff theorem (cf. Subsect. 3.4) predict the 
asymptotic values for SL and TL FF's to be identical at high energies, the experimental investiga- 
tion of this issue is confronted by the difficulties of disentangling the electric and magnetic FF in 
the TL region. If one assumes that the effective form factor could be an approximation of the TL 
magnetic form factor (cf. Eq. ([2])), one finds that it is larger than the corresponding SL quantities 
by about a factor of two. In order to perform this comparison, one has to assume that the TL 
magnetic FF is positive in the asymptotic region (cf. Subsect. 3.7 for a more detailed discussion). 



In Figs. (|9|and ([s]) the results of BaBar [9J and PS170 pjOj on the differential cross section da /dVl 
had already been presented. We remind that according to Eq. |4]the ratio \G'''^/G\^\ becomes accessible 
from the measurement of da/dVL. Fig. (15) shows the results obtained by both collaborations on this 
ratio. The experiment PS170 has measured \G^^/G\j\ in five energy bins below 2.1 GeV; BaBar has 
measured the ratio in six energy bins below 3 GeV. While the spectrum of the PS170 experiment seems 
to be compatible with the assumption IG^n/G^jl = 1, the BaBar spectrum shows a relatively large 
deviation from 1 for intermediate energies. 
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Figure 15: Only the PS170 (TU] and the BaBarjn] experiments have been able to measure the ratio of 
the proton's electric to the magnetic form factor. The results are incompatible with each other. A fit 
to the BaBar data is shown as well, cf. Ref. [9j. 



2.2.2 TL neutron electromagnetic form factor 

As mentioned above, only FENICE has measured so far the FF for the reaction e^e^ — > nn. In addition 
to the usual choice of using an effective FF (cf. Eq. (|2])), the FENICE collaboration has also analyzed 
the neutron FF by assuming \G%\ = in the one-photon cross-section, namely identifying the ratio 



3s 



with the magnetic FF of the neutron (see also Subsect. 3.2). This was motivated by the angular 
distributions measured in the experiment, which were compatible with a (1 + cos^^) shape, only. The 
data, hence, did not allow for an additional term proportional to sin^^, which would have indicated 
a significant contribution from the electric FF, see Eq. (g. In Fig. as in Ref. [131, the data are 
shown in terms of the previous ratio. 

Although the process e^e~ — )■ nn is measured with very low statistics, the FF is found to be 
systematically higher than in the case of the proton. This effect is not understood and needs an 
experimental confirmation, especially since naive quark models predict the TL proton FF to be higher 
than the neutron one. 



2.2.3 TL hyperon electromagnetic form factors 



The BaBar measurements of hyperon effective FF's are presented in Fig. 17 The results for the hyperon 
final states AA, AE*^ and E°E° are shown together with the proton FF. Although, the hyperon FF's are 
having significantly larger uncertainties confronted to the proton case, there is clear evidence that all 
measured FF's are featuring a strong threshold enhancement. Again, this behavior is not understood. 
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Figure 16: Measurement by the FENICE collaboration [13] of the energy dependence of the magnetic 
FF of the neutron In this analysis it has been assumed = 0, which was motivated by the 

angular distributions of the neutron-antineutron events. 



A nonzero relative phase between the electric and magnetic A FF's and Gj^^ leads to a 
polarization of the outgoing baryons with the polarization being proportional to sin0. The A polarization 
can be measured in the decay A — )■ pvr" using the correlation between the direction of the A polarization 
vector and the direction of the proton from the A decay. Due to the low statistics collected at BaBar, 
only a weak limit could be extracted for this phase: —0.76 < sin0 < 0.98. In Ref. |2H] also the 
|^e/^m| ratio was investigated in two mass intervals. Again, the low statistics did only allow for 
measurements with large uncertainties. IC^/G^^I is found to be 1.73!!^q57 mass range 2.23 GeV 

< Maa • c2 < 2.40 GeV and O.Tl+^J;?^ for 2.4 GeV < Maa ■ < 2.8 GeV, where Maa is the invariant 
mass of the AA system. 

2.3 Experimental outlook 

At present, the knowledge of the em FF's of the nucleon in the TL regime remains largely not understood. 
This statement appears obvious, given the various phenomenological puzzles mentioned above. The 
relevance of a better experimental data base for an improved theory prediction will be discussed in 
detail in Sect. [3j While experiments have been able to measure the shape of a so-called effective FF of 
the proton, |Ggff|, with a precision of a few percent at threshold, the uncertainties are much larger at 
higher momentum transfer. The neutron FF remains almost terra incognita. To receive a significant 
progress in our understanding of TL nucleon FF's, the future experimental programme, in our opinion, 
should concentrate on the following aspects: 

• Consolidate the existing FF measurements of the proton, and obtain first precision data for the 
FF of the neutron (as well as for hyperons); 

• Obtain statistically significant results for the ratio \G^/G^j\] 
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Figure 17: The energy dependence of the em FF of several hyperon states (AA, AE'' and S''E°) measured 
by BaBar [2S] is compared to BaBar's own measurement of the proton form factor [H]. 



• Measure the relative phase between the electric and magnetic Sachs FF's. 
Fortunately, in a world-wide effort such an experimental programme is currently pursued. All experi- 



mental techniques presented in Subsect. |2.1.1| will be used: the e~^e~ energy scan, the ISR technique, 
and the pp annihilation method. 

Very recently the VEPP-2000 electron-positron accelerator in Novosibirsk came into operation with 
cm. energies of up to 2.0 GeV and with a design luminosity of 10^^ cm~^s~^ at the highest beam 
energies. First preliminary results by the CMD-3 and SND collaborations concerning TL nucleon FF's 
have already been presented at conferences [53], both for the pp as well as for the nn final states. These 
first results indicate that - once the design luminosity will be achieved - we can expect approximately 
10^ pp and nn events produced close to threshold and the ratio \G^/G^\ to be measured with 10% 
accuracy. 

Furthermore, the second generation B- factory projects in Japan and Italy [51], which will start 
operation in several years from now, foresee an increase of luminosity of up to two orders of magnitude 
compared to the first generation. The ISR method will be applied at those facilities and the statistical 
errors for future FF measurements will be reduced by up to an order of magnitude. 

In the following, we will present in more detail two facilities for which a significant progress in the 
field can be expected in the near-term future, especially at intermediate and large energies, where an 
extraction of the \G^/Gfj\ ratio is of major importance. These are the BES-III experiment at Beijing 
and the PANDA experiment at the future FAIR facility in Darmstadt. 

2.3.1 BES-III at the e+e collider BEPC-II 

The BEPC-II e^e~ collider in Beijing is a next-generation r-charm factory with cm. energies corre- 
sponding to the mass of the charmonium resonances. Data taking started in 2009 with more than 
4 fb~^ of integrated luminosity having been collected since then. The world's largest statistics of J/ip, 
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Figure 18: The expected yield of pp events (left plot) and AA (right plot) using the ISR technique at 
a/s = 3.77 GeV at BES-III and assuming an integrated luminosity of 10 fb^^. The number of produced 
events for the BaBar kinematics, assuming an integrated luminosity of 232 fb~^, are also shown, for the 
sake of comparison. The PHOKHARA generator [l2l |13] has been used for these simulations. 



'0(3770) and ?/'(4040) decays is already available at BES-III, see Ref. [SS] for a detailed summary of 
the physics programme. The accelerator is designed to have the highest instantaneous luminosity on 
the '0(3770) resonance. Peak luminosities of 6 ■ 10^^cm~^s~^ have been achieved. This corresponds to 
~ 60% of the design value. For lower and higher energies, the instantaneous luminosity drops by up 
to an order of magnitude. Since the BEPC-II accelerator (differently from the B-factories) is designed 
to operate not only at one single beam energy, an excellent chance for a new precision measurement of 
the TL nucleon FF's over a wide energy range is given. 

Such an energy scan between 2.0 to 4.5 GeV is indeed scheduled for the coming years. In a first 
campaign, the step sizes will be approximately 100 MeV wide with a finer binning eventually at 2.25 
GeV and 3.0 GeV, where steps in the spectrum of the proton FF have been found (cf. Fig. (14)). 
Based on the current performance of the accelerator, the integrated luminosities per scan point can 
be estimated to be between ~ 10 pb~^ (at a/s = 2.0 GeV) and up to 100 pb~^ (below and above the 
■0(3770) resonance). Compared to the scan performed by the predecessor experiment BES-II (Ref. [T7]). 
this implies not only a significant increase of the number of scan points as well as of a broadening of 
the energy range, but also an increase in statistics of about two orders of magnitude. According to the 
current plannings, ~ 2, 300 detected pp events at y/s = 2.0 GeV and ~ 500 events at y/s = 3.0 GeV 
will be detected. Also concerning the systematic uncertainties, which used to be below 10% in the case 
of BES-II, a major progress can be expected, given the improved detector performance. The statistical 
error will however still be dominating the overall precision in most of the energy range. 

A determination of the em FF of the neutron seems feasible as well at BES-III. A recent measurement 
of the branching fraction B{J/ip — )■ nn) [56j with < 10% uncertainty has demonstrated the sensitivity 
of the BES-III detector for neutron and antineutron detection, which is the major issue in this analysis. 

The high integrated luminosities available at BES-III, will also allow for a programme of ISR mea- 
surements of TL nucleon FF's. Notice that the threshold region for NN production will only be 
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Figure 19: A sketch of the PANDA detector. 



accessible via ISR, given the fact that for technical reasons an operation of the collider below 2 GeV 
is most likely impossible. Background conditions appear to be optimal for runs taken on the ^^(3770) 
resonance, on which 2.9 fb"^ of integrated luminosity have already been taken and for which 10 fb~^ 
are expected for the coming years. A feasibility study ^7] shows, that the expected statistics at BES- 
III will be competitive with the existing BaBar measurement. The fact that the cm. energy of the 
BEPC-II collider is very close to the hadronic mass range of interest, leads to a radiator function Hs{s') 
(see Eq. ([s])) which is advantageous with respect to the measurements performed at BaBar. This effect 
cannot be overcompensated by the higher integrated luminosities available at the B-factories. Moreover, 
the fact that the ISR photon at BES-III is much less energetic compared to BaBar, opens the possi- 
bility for so-called untagged measurements, in which the detection of the ISR photon is not explicitly 
required. Since the differential cross section for ISR events increases significantly for low-angle photons, 
the untagged approach offers a very high statistics at a good signal-to-background ratio 

The expected BES-III statistics in terms of produced events using the ISR technique can be seen in 
Fig. (18) for the pp channel (left plot) and for the AA channel (right plot). An integrated luminosity 
of 10 fb~^ has been assumed in these simulations. For the sake of comparison, the event yield for the 
BaBar case for an integrated luminosity of 232 fb~^ is shown as well. Since the geometrical acceptances 
and detection efficiencies might be different between BES-III and BaBar, we only compare the number 
of produced event^ In order to facilitate the detection of ISR-produced nn events, a dedicated tagging 
detector [59] has been installed in the BEPC-II beam line at very small polar angles. This tagging 
device is intended to detect ISR-photons at very small polar angles. The precision of the future nn 
measurement remains however to be investigated. 



^The feasibility of untagged measurements was already proven in ISR-analyses of the TL pion FF at the KLOE 
experiment at DA<I>NE [58] . 

We remind that the fully available BaBar statistics is about a factor two larger than the one used in the published 
results of Refs. [21 [5S] and the one assumed in Fig. (18). 
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2.3.2 PANDA at the HESR antiproton beam at FAIR 



The HESR storage ring at FAIR [60j, which will start operation in 2018, will provide a high-intensity 
and high-resolution antiproton beam in the multi-GeV energy range. The beam will be scattered on a 
proton target (pellet or gas jet) within the PANDA detector, yielding luminosities of 10'^^cm~^s~^ and 
higher. This is an improvement of an order of magnitude compared to the luminosity, which used to 
be available at FNAL. 



The PANDA detector, which is shown in Fig. (19), consists of a barrel as well as a forward detector. 
The barrel detector resembles the geometry of typical collider experiments (see e.g. Fig. (10)) with a 
solenoidal magnetic field. The forward spectrometer reflects the fixed target concept of the experiment 
and contains a dipole magnet to allow for a momentum measurement of particles scattered at small 
polar angles. The PbW04 electromagnetic calorimeter of PANDA, the straw tube tracking system 
together with a dedicated DIRC Cerenkov detector setup will be of utmost relevance for a precision 
measurement of the pp — )■ e~^e~ process. With antiproton beam momenta between 1.5 GeV/c and 15 
GeV/c available at the HESR, the em FF of the proton in the range between 5 (GeV/c)^ and 28 
(GeV/c)^ can be accessed pTj . 

Detailed feasibility studies have been performed to prepare the measurement of the proton FF, cf. 
Ref. [52] • A serious background remains from two-body processes like pp — )■ tt+tt^ and pp — ?• vr^vr^, which 
happen to have up to six orders of magnitude higher cross sections compared to signal events. A special 
effort has indeed been made to provide realistic event generators for those processes (cf. Ref. [62]). In the 
pp — 7r°7r° background process, Dalitz decays of neutral pions can be rejected by appropriate angular 
cuts. The 7r~^n~ background channel requires very good PID capabilities and a good understanding of 
their experimental efficiencies. Kinematic fits help to reject the hadronic two-body background channels 
as well as other background processes. Taking into account all these analysis items, and by performing 
realistic Monte-Carlo based detector simulations, it has been found that the process pp — )■ e~^e~ can 
be measured over the full kinematically accessible range with a reconstruction efficiency of > 20%. 
The remaining background-to-signal ratios are expected to be as low as 0.1%. The measurement of the 
e~^e~ yield will be normalized to the integrated luminosity, for which a dedicated monitor system will 
be used at PANDA. 



In Figs. (20) and (21 ) the results of the above-mentioned feasibility studies are summarized. Fig. (20) 
shows the expected accuracy for the effective FF in comparison to already existing data. In the case 
of PANDA, an integrated luminosity of 2 fb~^ collected at 10 different values (i.e. scan points) 
was assumed; only the statistical error is shown. According to Fig. (20) we can expect a remarkable 
improvement with respect to the so-far most precise measurement from BaBar. Fig. (21) shows the 
accuracy (yellow error band), which is expected at PANDA for the ratio R = \G^^/G\.j\ (for this 
simulation it was assumed = |G^|). In the range between 5 (GeV/c)^ and 7 (GeV/c)^ the 
statistical error of the extracted ratio will be on the level of a few percent only. Above 14 (GeV/c)^ 
the extraction of the ratio is not feasible anymore due to the kinematic suppression of (cf. Eq. (|4])). 
It is obvious from Fig. (21) that a precision measurement of the |G^/G^^| ratio at PANDA will be very 
sensitive to distinguish between different hadronic models. Three of those theoretical predictions are 
shown in Fig. (21). The work corresponding to the red dashed line, cf. Ref. [65] (Brodsky, and Farrar), 
is predicting R = 1; the green solid line and the blue dot-dashed line are taken from Refs. [66] (lachello, 
Jackson, and Lande) and [67] (Lomon), respectively. Those models will be further discussed in Sect. |3] 

The PANDA collaboration has also suggested to measure the TL em FF of the proton in the 
unphysical region (g^ < 4Mp) by using the reaction pp — )■ e'^e'Tr", cf. Refs. [S3]. Furthermore, it 
has been demonstrated that TL measurements provide valuable information to search for two-photon 
effects, cf. Ref. [M]. 

As mentioned above and as worked out in Ref. [68], future directions in experiment should concen- 
trate on an experimental programme, which is aiming for a measurement of the relative phase between 
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Figure 20: Expected PANDA results on the TL 
magnetic FF of the proton and a comparison with 
effective FF's measured by various experiments 

(Refs. [la nn El [la EH M Ea Ea EI]) in the 

PANDA simulation it has been assumed: \G^^\ = 
IG^^I = |Ggff|. Figure reprinted from Ref. [22] with 
kind permission of The European Physical Journal 
(EPJ). 



Figure 21: Expected accuracy (yellow uncertainty 
band) at PANDA for the ratio of the magnetic 
to the electric FF. The results of the PS170 
|10j (squares) and BaBar [9j (triangles) collabo- 
rations are shown, as well as theoretical predic- 
tions, see text. Figure reprinted from Ref. [22] 
with kind permission of The European Physical 
Journal (EPJ). 



the electric and the magnetic FF's. In order to achieve this goal, a polarized target, or - as an alterna- 
tive - a polarized antiproton beam is requirecQ At PANDA, feasibility studies for a polarized hydrogen 
gas target are underway [7T] , while within the PAX collaboration [72] at FAIR solutions for antiproton 
polarization are being investigated. Spin filtering experiments are carried out in this context - so far 
with protons - indicating that this method seems to be a promising method (cf. Ref. [73]) for nucleon 
polarization. Hopefully this approach can be tested soon also with antiprotons. Double polarization 
experiments using both a polarized target as well as a polarized antiproton beam will be giving further 
insights concerning the relative phases between the magnetic and the electric FF's, cf. Refs. [UUl [70] 
for further details. In Sect. 3.2 we will be discussing in detail the opportunities which are given by the 
polarization degrees of freedom. 

The various plannings and future directions outlined in this Subsect. make us confident, that the ex- 
perimental progress expected for the coming years will bring us new and enlightening insights towards 
a quantum field theoretical understanding of the physics of TL nucleon FFs. We stress once more that 
future polarization experiments are a must. 



principle this phase can be also measured in the case of e^e^ — > pp annihilation using an unpolarized beam, 
but measuring the polarization of the outgoing proton. No experiment at an e^^e^ collider has however taken this into 
account. A polarimeter (e.g. made of carbon) would be required for this task. 
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3 Theoretical background 



The theoretical investigation of the nucleon em FF's in the TL region cannot be kept separate from 
the corresponding one in the SL region, given the complementary nature of the information on the 
nucleon internal structure, that can be accessed through the electron-nucleon elastic scattering and the 
annihilation processes NN -H- e^e~ or e"'"e~ -f-?- NN, all related via the crossing symmetry. Therefore, 
it appears useful to briefly illustrate some generalities on the nucleon FF's in the SL region. We remind 
that in a SL scattering process a virtual photon with a four-momentum = {u, q} hits a nucleon 
with four-momentum = {E = Mj^ + |pP, p}. The following relation holds for SL processes: 



= a;^- |q|^ < Q. 
3.1 Nucleon Form Factors in the SL region 

As it is well-known, in the SL region the nucleon em FF's can be investigated through the elastic 
reaction 

e- + N ^e'- + N' 



pictorially illustrated by the diagram in Fig. (22), where the one-photon-exchange approximation is 



shown. Within such an approximation, a virtual photon (g^ < 0) is exchanged between the incoming 
electron and the nucleon target, in a t-channel process. It is compelling to immediately mention the rel- 
evant role, which the polarization degrees of freedom played during the last decade in the experimental 
studies of the proton FF's (see Ref. |2] for recent reviews). In particular, by applying a novel experimen- 
tal technique, the so-called polarization-transfer technique, it has been discovered a totally unexpected 
feature of the electric FF of the proton: a fall-off faster than the dipole one, for —q^ > 2 {GeV/c)'^ (see 
below). 




Figure 22: The one-photon-exchange diagram for the electron-nucleon elastic scattering. The black 
bubble indicates the non-elementary nature of the photon-nucleon vertex. 

For < 0, the nucleon em current operator is a Hermitian operator, that behaves like a four-vector 
under Lorentz transformations, and, moreover, fulfills the constraints imposed by parity, time-reversal 
and current conservation. These general properties straightforwardly lead to the following covariant 
expression for the matrix elements in the SL region (see e.g. [Z1],[Z5] and [76]). 
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where i) Us{p) and Us'{p') are the initial and final nucleon spinors (normalized as Us'{p)us{p) = (5s', s), 
and ii) and are the Dirac and Pauli form factors, respectively. They are scalar functions that 
depend upon one scalar quantity, = 2(M^ — p ■ p'), namely the only, non trivial kinematical scalar 
constructed through p and p'. The second line in Eq. ^ has been obtained by applying the Gordon 
decomposition. It is worth noting that, in the first line of Eq. (|6]), the Dirac structure allows a physically 
grounded classification of the two terms. As a matter of fact, one has a helicity conserving contribution, 
that fulfills [7^7^^, 75] = and a helicity non- conserving term,that fulfills [7°cr'^'^, 75] 7^ 0, (let us remind 
that (1 ± 75)72 is the helicity projector for an ultra-relativistic case). 

The Dirac and the Pauli form factors take into account the non-elementary nature of the nucleon, 
and they are normalized according to the proton and neutron charges for and to the anomalous 
magnetic moments for , viz 

Ff (0) = 1 Fi"(0) = 

F^O) = 1.792847351 ± 0.000000028 ^2^(0) = -1.9130427 ± 0.0000005 ^ ' 

Notably, in the SL region, the Dirac and Pauli FF's are real functions. This feature can be deduced 
by using the Hermiticity of both and the Dirac structure (707^ and '~fo<^^'^)- One could infer the 
Hermiticity of J^^, by observing that the first-order approximation of the T-matrix is a Hermitian 
operator in the SL region, since in the t-channel no real particles can be present in the intermediate 
states. In particular, the unitarity condition for the first-order T-matrix, T^^\ reads 



I J^i^nr 5\p^ - q) [t«(^, q')] * T,«(^, q') = (8) 



As mentioned above, the vanishing value of the rhs is due to the unphysical nature of the intermediate 
states, \a). Then, the Hermiticity of follows from the Hermiticity of T^^^ oc Jem ■ je/q^- It is 
important to anticipate that, in the TL region, the unitarity constraint has to be implemented by 
taking into account the opening of real-particle production channels, since a s-channel reaction occurs. 
Therefore, the TL nucleon FF's are expected to be complex functions. 

Instead of using Dirac and Pauli FF's, it is possible to introduce proper combinations of them [TTflTS] 
as follows 

G^w') = F^m - TF^m , guq') = F^Q') + F^m , (9) 



where = _g2 > g, r = Q'^/AMl. In Eq. (|9|), G^{Q'^) and G^{Q'^) are the so-called Sachs electric 
and magnetic form factors, respectively. The physical motivation of considering these combinations 
is given by the analysis of the current operator performed in the Breit frame, where = |qp and 
p = — q/2 = — p'. In this frame, the nucleon FF's can be expressed in terms of three-dimensional 
Fourier transforms of functions to be interpreted as electric and magnetic densities, like in the non 
relativistic framework, viz 

pWM) = / ^'-"■-^Gf.M.dql'^). (10) 

where E{\q\) = ^/M^ + |qP/4. If -Edql) — 7- Mat, the non-relativistic limit is immediately recovered. 
Therefore, one could naively argue that the SL nucleon FF's yield information on intrinsic properties 
of the nucleon. Unfortunately, such a correlation should be carefully reconsidered, taking properly into 
account the boosts to be applied to the states of a composite system, like the nucleon (see, e.g., [75|[70]). 
Notice that, in the Breit frame, both initial and final states have to be boosted when the matrix elements 
of the current are calculated. 

Recently, a new physical interpretation of the nucleon FF's has emerged (see, e.g. |80j)- It is based 
on the field theoretical description of the constituents, inside the nucleon, given in terms of the so-called 
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Generalized Parton Distributions (GPD's) (see, e.g., [HI] for an extended review of the topic). Such 
a formahsm, devefoped within the Deep Inelastic Virtual Compton framework, leads to a description 
that fully satisfies the general principles, allowing to correctly treat the relativistic-boost effects, as 
well. It should be pointed out that GPD's depend upon three scalar invariants: x = k ■ n/{[p + p') ■ n], 
= {p' — p) ■ n/Kp + p') ■ n] and = {p' — p)"^, with i) k the average four momentum of the active 
constituent, before and after interacting with an external probe, and ii) n a light-like four-vector. 
Furthermore, it turns out that, in the Deep Inelastic Scattering (DIS) limit where the variable A plays 
an analogous role of q, one can reobtain the nucleon FF's from the unpolarized quark GPD's, putting 
^ = and integrating on x. To complete the analysis, two-dimensional Fourier transforms of the relevant 
GPD's in the Drell-Yan frame (a particular Breit frame where the + Qz = 0) are introduced. In this 
frame, the variable A'^ has only transverse components, i.e. A_|_ = q_|_, (notice that in the Drell-Yan 
frame = |q±P) and it is conjugated to the impact parameter b^. In conclusion, the main outcome 
is given by the determination of well-defined density distributions of the constituents in the transverse 
plane {bx,by}. For instance, the Dirac FF yields 



pf (|bj) = / ^ e^^-^- / dxH'^{x,0,A, 



1 



^ e^^-- F,^(g^ = |qj^) (11) 

where H^{x,0,A±) is one of the unpolarized GPD's. Summarizing, can provide information on the 
density distributions in a transverse plane, once the effect of the boosts on the initial and final states of 
a composite system are correctly taken into account. This can be implemented by using the light-cone 
description of both the operator current and the nucleon states (see, e.g., [SI]). 

The last issue to be mentioned, that will be relevant also in the TL region when a model builder has 
to constrain his own approach, is the asymptotic behavior of the nucleon FF's, as dictated by QCD. By 
considering both dimensional scaling laws and helicity conservation, perturbative QCD (pQCD) yields 
the following asymptotic behavior for the nucleon FF's (see, e.g., [52]) 



1 



In 



1 —1 



A2 



1,2 (12) 



where 7 is the anomalous dimension (7 = 2.148, 2.160, 2.173, depending upon the active number of 
flavors Nf = 3,4,5). It should be pointed that the asymptotic behavior of the FF's is dictated by the 
perturbative gluon exchange, that provides the mechanism for sharing the momentum transfer among 
the constituents. In order to accomplish this, at least two massless gluons must be exchanged, and each 
gluon propagator contributes with a factor l/Q^. As to , the spin-flip effect adds an extra factor 

To complete this snapshot of the SL region, it is impossible to avoid few words on the experimental 
technique that has given a new, strong impetus to the investigation of the proton FF's in this region. 
In particular, a technique based on the measurement of the polarization of the final state proton has 
experimentally proven that the proton charge FF strikingly decreases faster than the dipole FF [2] . It 
is also worth noting that in the TL region, since the nucleon FF's become complex functions, a NN 
pair, produced in the final state, naturally acquires polarization degrees of freedom, even in absence of 



polarized leptons in the initial state (see Subsect. 3.2 and |68]). 



In the SL region, when a polarized electron hits an unpolarized nucleon target, some polarization 
is transferred to the recoiling nucleon. Such a polarization can be measured with a suitable secondary 
analyzing reaction. The polarization-transfer method (see [2] for more details) is based on the mea- 
surement of the transverse and longitudinal components of the recoil proton polarization, in the elastic 
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reaction 

e-+p^e'~+p' , (13) 

where the incident electron has a longitudinal polarization ±/i along the direction of the momentum 
transfer, q = z (indeed one could consider reactions with more complicate polarization set-up [83])). To 
understand why the final state proton gets a polarization, one should first recall that the polarization 
of a fermion (with four-momentum p and spin s) with respect to a given direction, n, is obtained from 
the following expectation value: 

('^ ■ n) = --, \-, ^u\p,s) i ^ rr^ i,] ^(P'^) = —( TTT V u{p,s)-i^ |4m(p,s) = 

w{p,s)u[p,s) \ U (T ■ n J u[p, s)'j^u[p, s) 

^777 7 75 ?i = —pr s ■ n, (14) 



E V 2M/V 2 W E 



where = {0,n}. In order to determine the polarization of the final proton, Eq. (14) has to be 
properly weighted by using the polarization density of the final state, i.e. by exploiting the transition 
probability to scatter the final proton along a given direction. To this end, one has to use the density 
matrix of the final nucleon, formally written as p/ = AipiM\ where the matrix Ai is the T-matrix 
relevant for the process under scrutin}|^ (see, e.g., [83J for a discussion in one-photon approximation 
within a relativistic framework). It is worth noting that the density matrix of the final nucleon can lead 
to a polarization in different cases: i) in case the initial nucleon has a polarization, ii) if the incoming 
beam has a polarization, or iii) a complex-valued FF's of the nucleon are present in the current (see [HI] 
for the non relativistic case, and [68] for the relativistic one). In conclusion, the above mentioned 



weighting of Eq. ( 14 ) amounts to obtain the polarization of the final proton by evaluating the polarized 
differential cross section, where an average on the spin of the initial proton is performed (this means an 
initial density matrix purely scalar), and the polarization of the incoming beam is taken into account. 
Then, the experimental polarization, P = {Px,Py,Pz}, of the recoiling proton is given by (see, e.g. [SS] 
and [Hi]) 

P-h^ = Kj L,^W^: (15) 

where da/dfl is the unpolarized cross section (representing the normalization of the polarization den- 
sity matrix in the final state), Kf a proper factor, L^,^ the leptonic tensor with both symmetric and 



antisymmetric terms (since the electron beam is polarized in (13)), and the proton tensor is given by 



= Tr \tl±^ J^tl±^ r 

^ ^ 2Mp 2 ^ 2Mp P 



It turns out that the proton in the final state of the reaction (13) gains a polarization in the scattering 
plane, {x,z), taking q = z. The measured components of the polarization vector allows one to extract 
the ratio of the proton FF's as follows 

Gl P^Ei + Ef f9e\ 

-^ = -— ^ tan — (16) 

Concluding this Subsection, for the sake of illustration, an overview of the status of the analysis of the 



nucleon FF's in the SL region is presented in Fig. (23). It is also shown a recent theoretical calculation 



[7], based on a relativistic constituent quark (CQ) approach, that has been extended to the TL region. 



as well (see Subsect,3.9 ). The interested reader can be referred to, e.g., Ref. [2] for a comprehensive 
and updated discussion of the experimental measurements. 

^In a non relativistic framework, one parametrizes M = g — ih ■ cr, with g and h complex functions, and Pf = 
(da/dn) [l + P-cr]/2. 
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Figure 23: Space-like nucleoli form factors vs. 



-q^. Left panel: the proton form factors. 



Right panel: the neutron form factors. Data from the compilations in [2]. The dipole form factor 
GD{\(f\) = [1 + k^l/lO.T^GeV/c)^)]"^ yields the reference line of the FF's, when indicated. The solid 
lines are the calculations, with four free parameters, performed within the relativistic Constituent Quark 
Model of Ref. [7]. (Adapted from Ref. [7]). 



(17) 



3.2 Nucleon form factors in the TL region 

The reaction where a NN pair annihilates into charged leptons, viz 

N + N + r 

and the time-reversed process, viz 

t + r ^N + N (18) 
represent a direct source of information on the nucleon FF's in the TL region. For the sake of con- 



creteness in what follows it will be considered the e^e annihilation. In Fig. (24), the one-photon 



approximation of the reaction ( 18 ) is depicted. 



In the hadronic annihilation process, a virtual photon with invariant square mass > can mate- 
rialize into a quark-antiquark pair, since it can become active the hadronic part of the virtual-photon 
state, shortly given by 



I7*) = \point — like) + |e^e ) + |/i 11^) + .... + \qq) + 



+ \qq99) + 



leptonic term 



hadronic term 



(19) 
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Figure 24: The one-photon-exchange diagram for the electron-positron annihilation into a NN pair. 
Notice that a physical NN pair is produced only for > 4M^. As in the SL diagram, shown in 
Fig. (22), the black bubble indicates the non-elementary nature of the photon- nucleon vertex. 



An honored, and still very alive, approximation of the photon vertex is the vector dominance model 
(VMD) (see, e.g., [3J for a recent review), that describes the hadronic part of the photon wave-function 
through a direct coupling between the virtual photon and a vector meson (with the proper mass, spin 
and parity). Such a model certainly allows one to construct a very effective scheme of approximation, 
that has to be completed by a description of the hadronic decay of the vector meson, if one is interested 
in the hadronic FF's (i.e. 7* — VM — )■ hh). In view of the quark content of the hadronic term, one 
should also consider a direct coupling to a qq pair, and therefore more complicate production channels 
(as it will discussed in what follows). With these first physical insights in mind, we can start to deal 



with the analysis of the matrix elements of the nucleon current for describing reaction (18). 



The investigation of the process (18) (as well as (17)) takes advantage of the crossing symmetry, and 
therefore one can exploit the previous analysis of the elastic scattering in the SL region. In particular. 



the matrix elements of the nucleon current operator involved in the reactions ( 18 ) are written as follows 
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(20) 



where = {p' + p)^ = AE'^, = {E,p} and p''^ = {E,—p}, in the center of mass. In the TL 
region, the current operator is involved in the transition from the vacuum to a state with a nucleon- 
antinucleon pair, that becomes real beyond the proper threshold, given by = 2M]\f. This means 
that production channels are now opened, and therefore the current operator is not more a Hermitian 
one. This can be verified from the unitarity condition, Eq. ([s]), applied to this case, namely to an s- 
channel process. As a consequence, as it can be also recognized from a perturbative analysis of the FF's 
themselves, the nucleon FF's become complex functions, with a characteristic multi-cut structure due 
to the free propagation of the relevant hadronic pair, beyond the proper production threshold. Indeed, 
if we expand the analysis to include the whole hadronic production, one can recognize that different 
channels can be opened, with different energy thresholds. In the interval < < 4m^, there is no 
purely hadron production at all, while for increasing values of g^, up to g^ < 4M^, one meets channels 
that contributes to the production of a virtual NN pair. One has in order : i) the isovector two-pion 
channel, that starts at g^ > 4m^, where m,r is the pion mass, ii) the isoscalar three-pion channel, that 
begins at g^ > 9m^, iii) the (p meson production, beyond 1 GeV^, etc. The physical threshold for the 
NN is g > 4M^ and the region 4m^ < g^ < 4M^ is called the unphysical region, when we are dealing 
with the reaction (18): it contains many interesting information, particularly near the NN threshold 



see Subsect. 3.3). It should be pointed out that the opening of more production channels beyond the 



NN threshold generates new overlapping cuts in the FF's. As we will see below (Subsect. 3.4), one of 
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the main issue in the general treatment of the TL FF's is precisely the evaluation of the discontinuity 
of the FF's across the cuts, namely the imaginary part of the FF's. 




By adopting the one-photon approximation and performing the suitable traces, one gets (see also 
Eq. the following differential cross section for the A^iV-pair production, in the CM frame, [86] (for 
the cross section of the reaction pp — )■ e~^e~ see Ref. [87], paying attention to the flux factor different 
from 1, in this case) 

(^T =^ C,{q') [\Gim\l + cos'0) + \G^{q')\' Isin'e] (21) 



where a = e^/(47r) ~ 1/137, 6 is the angle between the direction of the incoming electron, taken as 
z-axis, and the produced nucleon (see Fig. ( 25|), (3 = a/1 — 4:Mff/q'^ = \p\/E the absolute value of 



the velocity of both N and (with E = \/Mfj + |p|^), r = q'^/iMj^ (notice the difference with the 
definition in the SL region, below Eq. (IqI)) and Cat is a constant defined as follows 




CN{y) = { (22) 



with 

Mpan 

y = 2 w- 

q /3 

(notice that the expression of the variable y contains a factor 2 at variance of what one can find in 
Ref. [9j where a mistyping is present [HB])- In the case of the pp-pair, Cp{y) is the s-wave Sommerfeld- 
Gamow factor that takes into account the QED leading-order correction to the wave-function of the 
charged pair, and results to be proportional to |\E'pp(0)p, where ^'pp(O) is the relative wave-function in 
the continuum. This factor specifically affects the cross section near the threshold, where the relative 3- 
momentum of the charged hadronic pair is small (see, e.g., [891 190]) and therefore /3 — )■ 0. In particular, 
it makes the cross section different from zero at threshold, since the factor (] in the denominator of 
Cp{y) ~ y balances a factor f3 that comes from the phase-space calculation. As a final remark, it is 
worth noting that the standard non relativistic expression of the variable y is 

y = = zair - — - , 

Vr IpI 
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where Vr is the relative velocity, |p| the relative momentum and firid the reduced mass (Mp/2). It can 
be heuristically generalized to the relativistic framework [9J as follows 



2E IpI 
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/3 



obtaining the expression of the variable y in Eq. ([22]) (notice that c = 1 in the above equations). 



In Eq. (21), as in the SL case, one has introduced the Sachs FF's, given by 
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(23) 



Notably, at the threshold, i.e. q^ = 4M|r, electric and magnetic Sachs FF's become equal by construc- 
tion. This issue has been experimentally investigated for the proton case, by two different collaborations. 



As discussed in Subsect. 2.2, the ratio |Gg(g^)/G^(g^)| has been extracted from the angular distribu- 
tions (cf. Eq. (21)) measured in a wide range of momentum transfer, 2Mp < q < 3 (GeV/c), by the 
BaBar collaboration p. For q up to 2.2 (GeV/c), the ratio results to be significantly greater than unity, 
in disagreement with previous data, obtained by the PS170 collaboration ^Q\, that are compatible with 
the expected threshold value ~ 1. 



The total cross section in one-photon approximation is obtained from Eq. (21) by integrating over 
the solid angle . For the sake of legibility, here we rewrite Eq. ([T]), but with some obvious changes. 
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Such an expression suggests to define an effective nucleon form factor in the TL region by dividing 



1/^2 



the actual total cross section shown in Eq. (24), by the point-like one, obtained from Eq. (24) putting 



1. Then, one can write 
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where 



The total cross section becomes 
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(26) 



(27) 



The point-like cross sections for proton and neutron reduce at threshold, i.e. r = 1, 

2 3 

<o.nt(4M|) = ^ = 0.848 nb <_,(4M„2) -> 

where for the neutron C„/3 — )■ for — )■ 4M^. In Fig. (26 ), for illustrating the effect of the Sommerfeld- 
Gamow factor and the behavior close to the threshold, a^j„j(g^) for both proton and neutron, is shown 
in the range of the present-day experiments. One should notice i) the very steep increasing between 
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Figure 26: Point-like nucleoli cross section, see Eq. (26) and text, vs q . Solid line: proton. Dashed 



line: neutron. Notice the effect of the Sommerfeld-Gamow factor, Eq.(22), at the threshold of the pp 
production. 



the thresholds, g^^^ ~ 3.52 (GeV/c)^ and q^^/^ ~ 3.53 (GeV/c)^ respectively, and the maxima, around 
q"^ ~ 4.9 (GeV/c)^, and ii) the effect of the Sommerfeld-Gamow factor that makes non vanishing the 
point-like cross section of the proton, at the threshold. 

Coming back to the definition of the effective FF, that in one-photon approximation is a combination 



of electric and magnetic Sachs FF's, it appears quite natural to compare Eq. (25) with the corresponding 
quantity defined through the experimental cross section o"g+g__j.jYjy, namely Eq. (j2]). It should be pointed 
out that in both numerator and denominator in the middle of Eq. ^ there are quantities experimentally 
measured, and therefore one has an operative definition for G^(g^), a part the model dependence due 
to the Sommerfeld-Gamow factor. Moreover, in analogy with the microscopic interpretation of the SL 



FF's as Fourier transform of electric and magnetic spatial distributions (cf Eq. (10) with the caveat of 



Eq. (11 )), one could develop a simple physical picture of the TL effective FF. Indeed, one could consider 



that it is microscopically related, through a Fourier transform in time, to the transition amplitude from 



a virtual-photon state (cf Eq. (19)) to a hadronic state. With this in mind, it becomes clear that in the 
TL region one can explore the time/energy structure of the hadronic Fock components of the photon 
state. At this point, a final short comment about (j'f^(g^) seems appropriate. One should notice that, 
unfortunately, in the literature one can find a shortcoming in the presentation of what the experiments 
provide. As a matter of fact, sometimes it is claimed that one measures the nucleon magnetic form 
factor assuming = but this is not the case, unless r = 1 (i.e. at the threshold). 

However, it should be clarified that the proton database in the TL region, is entirely constructed by 
introducing the effective proton form factors as defined in Eq. ([2]). For the neutron case, the Frascati 
data points [13] are shown i) by using Eq. ^ or ii) by assuming \G\{q^)\ = 0, that amounts to define 
for the neutron 



\Gfrr.ig')\ = j "'"'rX^''^ (29) 



47ra2/3 
'3q 



2 



In any case, in the experimental papers the procedure adopted for extracting the data is always properly 
indicated. 

In the TL region, since the nucleon FF's are complex functions, one has extra degrees of polariza- 
tion for the produced NN pair, even without a polarization of the incoming beams. This was already 



33 



explained above, in the brief introduction to the polarization-transfer method in the SL region. For in- 
stance, one finds that a produced proton jSHl EH EB] gains a polarization perpendicular to the scattering 
plane (recall that the chosen z-axis is directed along the incoming electron beam) and, in one-photon 
approximation, it is given by 



.m(2^) '^m{Gl{q')Gl,^{q')] 



D ^ 



D 

where ^e{m) is the phase of the complex- valued electric (magnetic) Sachs FF and D is defined as: 

D = \GlM')\\^^cosH)^\GW)f— 



(30) 



(31) 



The other two components of the polarization, and P^, lie on the scattering plane and are different 
from zero only if the electron beam has a non vanishing longitudinal polarization, Pe, viz 



- P. 



P. = P. 



2cos{9) \GMp{q 
D ^ 



D ^ 

2m2 



(32) 



As one can see from Eqs. (30), (31) and (21), information on absolute values and phases can be extracted 



by measuring both the angular distributions and the normal (to the scattering plane) polarization. 



Indeed, one could consider also different cases, e.g. a polarized antiproton beam in the reaction (17) [72] . 
in order to extract double-spin polarization observables, that allows to access '^t{G^G^^^ (see Refs. [69] 
and [To] for details). 

The previous expressions of both cross section and polarizations have been obtained by considering 
only the one-photon-exchange approximation. The formalism for including the two-photon-exchange 
effects can be found, e.g. in Ref. [92] (for a recent review of the effects of the two-photon exchange 
in the SL region see Ref. [93]), while Refs. [9ll |95] give a first quantitative investigation of the effect, 
showing a reassuring smallness (few percents) of the corrections. 

Finally the asymptotic behavior has to be the same as in the SL case (cf. Eq. (12)) since it is 
governed by pQCD (for an extended analysis of this issue see also Ref. [55]). 



3.3 The A^A^ production threshold 



The study of the reaction (18) near the threshold provides many relevant insights in the reaction 
mechanism that governs the transition from the unphysical to the physical regions. 

For the proton case, the PS170 collaboration [TD] observed for the first time a steep slope of the 
effective FF when g° slightly increases from the threshold value, and more recently the BaBar collab- 
oration [9] has accurately measured and confirmed the near-threshold peak. As to the neutron, the 
FENICE collaboration [T3] (see also [T2] for the proton measurements) made the first and (till now) only 
measurement close to the threshold, showing also in this case, but with large error bars, an increasing 
value for g° ~ q^^j^. 

The sizable and sharp rising of the cross section close to the pp production, given the experimental 
efforts devoted to its investigation, has driven a lot of theoretical studies, pointing to elucidate the 
origin of such a surprising result. It is useful to remind that for the pion, one meets a more soft 
behavior, since a huge peak related to the p-meson production, around 770 MeV, appears far beyond 
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the production threshold of a tt+tt^ pair, q^^i^ = 2m^. Even if the proposed mechanisms for explaining 
the proton case are different, all the approaches have the common aim to clarify and constrain the NN 
interaction. In particular, the steep slope has been explained i) by the pp final-state interaction acting 
near the threshold (see, e.g. [97] and references quoted therein), ii) by the tail of a pp bound state (see, 
e.g. |98] for a recent investigation of this issue by using the NN Paris potential and reference quoted 
therein) or iii) by a narrow meson resonance [L3\, just below the threshold. For instance, in Ref. [99] a 
complex scattering length in a ^5*1 pp state is obtained from the Jiilich NN potential |1UU] . allowing to 
achieve a good description of the data. In Ref. [97] , the field of investigation is enlarged to cover not 
only the proton case but also nn and other BB systems, by exploiting a phenomenological approach 
beyond the scattering-length approximation, pointing to a resonant behavior near the threshold. The 
pp bound state appears a very intriguing possibility, but such an explanation as mentioned above, is 
in competition with less exotic effects. However, if this is the case, one could address the issue via 
reactions like e~ + e"*" j + pp or even p + — e~ + e"*" + n as investigated in jlUl] , by taking into 
account the off-mass-shellness of the nucleons inside the target nucleus. 

The neutron threshold deserves some attention, in view of the existence of only one data set [13] . 
Those data posed a question about the applicability of a naive pQCD analysis to the qq production 
involved in the proton and neutron FF's. By using the FENICE data for both proton and neutron, one 
has |G^(4M^)|/|G'p'^(4M^)| ~ 1.4, while a naive perturbative description of the e~^e~ annihilation [lU2j 
leads to |G'^(4M^)/Ggg(4M^)| = \ed/eu\ = 1/2, namely almost a factor 3 of difference. Fortunately 
the new data by the BaBar collaboration [9] suggest a path that reconciles experiment and naive 
expectation, since by considering the BaBar proton data near threshold, the difference between the 
experimental ratio and the naive pQCD prediction shrinks the difference to a factor of 2, but with a 
large error band, due to the uncertainties of the FENICE neutron data. 

3.4 Dispersion Relations: an overview 

A truly microscopic calculations of the nucleon FF's have necessarily a prerequisite: i) a knowledge of 
the substructure of the nucleon and the interactions between the constituents, and ii) a non-perturbative 
framework where elaborating the calculations. Such an observation makes clear how much challenging 
be the task. Historically, coming back to the fifties, the study of the nucleon FF's had great benefits 
from the extension of the investigation to the whole complex plane of the relevant variables. In this 
way, as briefly illustrated below, the non-perturbative regime of the strong interactions, based at that 
time on nucleonic and mesonic degrees of freedom, could be taken into account, by exploiting both i) 
information gathered within a perturbative (diagrammatic) analysis of the needed amplitudes and ii) 
more fundamental properties, like causality, embedded and tightly related to the analytical behavior 
of the observables. As it is well-known, such an analytic method has its origin in the Kramers-Kronig 
relation in optics \103\ . 

Local properties of the field theories, as dictated by the microscopic causality, can be immediately 
translated into the analytic behavior of the scattering amplitudes (see, e.g. [76]), and in turn into the 
analytic behavior of the general Green functions of the interacting fields, if one adopts the Lehmann- 
Symanzik-Zimmermann (LSZ) reduction formula [1U4] for the S-matrix elements. 

The causality constraint imposes a finite value for the speed of any signal. This constraint can be 
rephrased into the request that space-time regions, separated by space-hke distances, cannot commu- 
nicate each other and therefore one has to require proper operatorial relations between fields, viz 

[(j){x'),(j)\x)] =0 {^{x'),^^{x)} =0 for (x'-x)' = c'(t'-t)2-(x'-x)2 < 0(33) 

where 0(x) and ip{x) are the local bosonic and fermionic fields of a given theory, respectively. In 
turn, the commutation properties establish well-defined analytic behavior of the Fourier transform of 
combinations of the field themselves, analogously to the basic investigation performed in classical optics 
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by Kramers and Kronig, who analyzed the propagation of a wave front and found the constraints 
induced by causahty on the diffractive index (i.e. the em response of the composite system). A causal 
behavior means that the response of the system will be vanishing for t' < to, where to is the instant 
when the system starts to be perturbed. If one takes into account such a constraint, one can rigorously 
deduce that the Fourier transform of the system response must be analytic in the upper-half plane 
(see [76] for details). Then, an integral relation between the real and the imaginary part of such a 
Fourier transform follows. For instance, if one has a complex function F{z) = ^eF{z) + i'^mF{z), 
analytic in the upper-half plane, one can write the following Cauchy integral 

F{z) = — [ dz' (34) 
^ ^ 27ri Jc z' - z ^ ^ 

with a contour, C, composed by the real axis and a semi-circle at infinity, in the upper half-plane. If 
F{z) properly behaves for \z\ — )■ oo, namely the contribution from the arc at infinity is vanishing, one 
gets a simple expression 

Fiq) = lim F(q + ie) = / dq ^-^ — (35) 

otherwise one has to add the contribution of the semi-circle. Moreover, recalling the formal expression 

1 V 



q' — q ^le q' — q 



± m6{q' — q) 



where V indicates the Cauchy principal value, the following relations, called dispersion relations (DR's 
can be straightforwardly deduced 

SeF(,) = r M W) (36) 

3mF(g) = - -V r dq' (37) 
7^ J -00 1 -Q 

Therefore, the knowledge of the whole function F{q) can be obtained from the knowledge of only one 
contribution: the real part or the imaginary one. The name dispersion relations can be traced back 
to the work of Kramers and Kronig |103j . where the diffractive index was investigated, and relations 
between the real part (governing the dispersive behavior of the light) and the imaginary one (controlling 
its absorption) were established, by using the causality principle, as mentioned above. 



If the assumed quick fall-off of F{z) does not hold, and therefore the integrals in Eqs. (36) and (37) 
are not well-defined, then one can write down new DR's, by introducing a proper subtraction term. For 
instance, for F{0) 7^ 0, one can speed-up the fall-off by subtracting F{0)/z and repeating the same 



steps that lead to Eqs. (36) and (37), but for 

Fiz) - F(0) 



obtaining, e.g., for the real part 



3fJeF(g) _ 3fteF(0) ^ 1^ QmFjq') 



Q q n J_^ q'iq'-q) 

since V dq' l/[q'{q' — q)] = 0. Such a subtraction procedure can be applied with the obvious changes, 
e.g. subtracting F'{0)/z from [F{z) — F{0)]/z'^, as many times as it needs in order to make vanishing 
the contribution at infinity, in the integrals. 



^These relations are also known as Hilbert transformations between dieF{q) and QmF{q). 
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In view of the application to the nucleon FF's, it is interesting to briefly discuss the case when F{z) 
has a cut for go < < oo, and it is real for — oo < ^ez < go- The Schwarz reflection principle, 
F{z*) = F*{z), applied along the open interval — oo < ^ez < go, allows one to analytically continue 
F{z) in the lower half-plane. Then a different path C can be considered for writing down the Cauchy 
relation (34). In particular, the integration path runs along the lower rim from infinity to go and 
along the upper one in the reversed direction, and it has to be completed by a circle at infinity (to be 
integrated in the counter-clockwise direction). Then, assuming for simplicity that the contribution to 
the Cauchy integral for |z| — )■ oo is vanishing (a subtraction procedure can be always implemented), 
one gets 

f (.) = ^ hm r = i r w) 

In conclusion, the function F{z) is determined for any z by the discontinuity across the cut. The 
relevance of this relation is given by the possibility to link the imaginary part of the nucleon FF's in 
the TL region to physical amplitudes, as discussed in the following Subsection. With simple steps, one 
can also write 

^eFiq) = lim ^eF(q + te) = -V / dq' ^ (40) 

^^0+ ' ' vr ^ ^ g'-g 

Such an equation is a generalization of Eq. (36), where the analyticity of F{z) was assumed only in the 
upper half-plane. In the case of a slow decreasing of F{z) for \z\ — )■ oo, one can repeat the subtraction 
procedure as in the previous Eq. (p8|), obtaining the following one-time subtracted expression 



Fiz) = F(0) + - lim r dq' ^^;^^^^'\ , (41) 

In closing this Subsection, it is important for what follows to mention one of the theorems that 
can be deduced from the Phragmen-Lindelof principle that allows to extend the maximum-modulus 
principle to unbounded domains |105] . The theorem, relevant for putting in relation the asymptotic 
value of the FF's in the TL region and in the SL one, states: 

// an analytic function F{z), with z = re'^'^, is such that i) F{z) a as z ^ oo along a straight line 
starting from the origin and F{z) —^basz^ — oo along another straight line, again starting from the 
origin, ii) it is regular^nd bounded in one of the two sectors between these two straight lines. Then 
F[z = oo) = F[z = — oo) = a and F{z) — )■ a uniformly in that sector as \z\ — )■ oo. It is important to 
note the request that the analytic function must be bounded in one of the two sectors. 

Summarizing the analysis in the complex plane, one should mention four items: i) causality de- 
termines the analytic behavior of an amplitude involved in a given physical process, ii) the Schwarz 
reflection principle allows to analytically continue a function from one side of the plane to the other side, 
if along a part of the real axis it does not contain singularities, iii) in presence of a cut, the function can 
be expressed in terms of the discontinuity across the cut, iv) a theorem from the Phragmen-Lindelof 
principle leads to relate the asymptotic behavior of the function for z — t- ±oo. One of the first discussions 
of the application of the Phragmen-Lindeloff theorem (previously used for proving the Pomeranchuk 
theorem and its generalizations) to the em nucleon FF's, putting in relation the asymptotic values in 
the SL and TL regions, can be found in Ref. 



3.5 Dispersion Relations and the Nucleon Form Factors 

Refs. [U 1106] (see also [76] , [H] and |107j for more general discussions) yielded a first description of 
how the previous general properties, particularly causahty, can be exploited for modeling the nucleon 
FF's in both the SL and the TL regions. Noteworthy, the same approach can be applied to any hadron. 

regular function is analytic and single- valued in a given region. 
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In the seminal work of Ref. [f4j, dispersion relations for the SL and TL nucleon FF's were determined, 
after analyzing the matrix elements of the em current operator (needed for the first order T-matrix) 
and using the LSZ reduction formahsm |104] for the S-matrix elements (see, e.g., [76]). A direct and 
fundamental consequence of such an investigation was the possibility of relating the imaginary part of 
the FF's to matrix elements relevant for other physical processes, that, in turn, can offer experimental 
inputs for an overall description of the FF's themselves. Let us briefly illustrate the issue in the TL 
case, by considering the following matrix elements of the em current operator 



^^'(2vr)3 {0\Jif^JOmp,sWip',s'))=vAp') 



Ml 



ZlVi]\[ 



Us{p) 



= ^^i^rrf' J d'x e-^^'- vAp'W\r+{SN{x)JUm^{P, s)) (42) 
where, due to the LSZ reduction, one has 

r+{5^(x)JL(o)} = e{t)Sr,{x)JU^) + e{-t)j^^{Q)SN{x) = 

= e{-t)[JU^),SN{x)] + SM{x)r,^{()) , (43) 

with Sn{x) the source of the nucleon interacting field, i.e. 

{if -m)^N{,x) = Sn{,x) (44) 

It should be recalled that Sn{x), like ifjN^x), lowers the baryon number by a unity. 

An important step for making more transparent the investigation of the analytic behavior of the 
FF's, is the substitution of the time-ordered sum 7+{S'Ar(x) J^^(O)} by the commutator [^^^(0); Siq{x)], 



that allows one to implement the causality (cf. also Eq. (33)). This can be accomplished by noticing 



that the matrix element (0|S'Ar(a;) J^„(0)|A/'(p, s)), generated by the rightmost term in the second line of 



Eq. (43), does not contribute to the matrix element of the current, Eq. (42). As a matter of fact, since 

(if — m){Q\'tp]^{x)\one nucleon) = 

the matrix element {0\SN{x)\one nucleon) must vanish. Therefore, no contribution can be generated 
by a one-nucleon intermediate state, when a physical, complete basis is used for spanning the matrix 
element {0\SN{x)J^^{0)\^^{p,s)), viz 

{Q\SAx)Ji^^iO)\Xip,s)) = -i^ Yl {0\SNix)\pi;I,B = l){prJ,B = l\J^^mM{p,s)) = 

Pl',I^nucleon 

= 7^ E (0|5;v(0)b,;/,i? = = 1| JL(0)|Ar(p, .)) (45) 

pj-Jj^nucleon 

where Ylpj-i Cleans an integration over the four-momentum and a sum over the discrete quantum 
numbers. It should be pointed out that the intermediate states spanning the matrix element in the Ihs 
of Eq. ( [45| must have baryon number B = I, since both the em current conserves the baryonic number 
and Sn{x) annihilates a nucleon in the initial state. Moreover, pj > Ml (the nucleon is the lowest 
B = 1 state), and therefore the intermediate states must contain multi-hadron systems with B = 1. 



Finally, since in Eq. (42) {p'Y = Ml, one immediately deduces that (0|S'Ar(a;) J^^(0)|A/'(p, s)) does not 
contribute to {0\J^^{0)\Af{p, s)J\f{p' , s')) (see also [HI] for details). As a matter of fact, one gets 

I d'x e-^^'-^ (0|^^(x)JL(0)|Ar(p,s)) = ^\P+Pi) 

pi\I^nucleon 

x(0|^^(0)|p,;/,i? = l)(p,;/,i? = l|JL(0)|Ar(p,s)) =0 (46) 
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Therefore, the matrix element in Eq. (42) can be rewritten as follows 



'EE' 



E 



M, 



(47) 



As mentioned above, the presence of the commutator makes it possible to exploit rather directly the 
causality, since [J^^(O), S'Ar(x), ], is vanishing for space-like separations, i.e. < (given the commu- 
tation rules of the local fields). Then, one has to deal with two theta functions in the four dimensional 
integral (47): 9{—x'^) and 0{—t). A general discussion, that involves a careful analysis of the dependence 
in exp{—ip'-x), where p'-x = i?'t — p'-x and |p'| = a/ E''^ — M^, can be found in [HI]. Here, following the 
simplified presentation of Ref . [4] , the case p' = is assumed, and the analytical dependence upon only 
E', taken as a complex variable, is considered. Given the function 0{—t), the analytic structure of the 

QmE'), and this leads to argue that the matrix elements 
^eE' + i '^mE'. Straightforwardly, one can 



integrand in Eq. (47) is governed by exp{—\t\ 
of the current is analytic in the upper half-plane of E' 
repeat the steps illustrated in the previous Subsection and construct dispersion relations for the nucleon 



EE's. Indeed, if one would use the dispersion relations (see Eq. (39) or its subtracted version, Eq. (41)) 



one has to fully determine where the singularities are located, and in this task the help comes from 
the analysis in perturbation theory, namely by considering the Eeynman diagrams contributing to the 
EE's. Such a study (see, e.g., [70]) shows that in the TL region one meets multiple cuts: a first branch 
point located at s = 4m^, i.e. in correspondence with the threshold of the two-pion production, that 
contributes to the isovector part of the nucleon EE's (in a VMD framework [3], one has the production 
of a virtual p meson, that decays in two pions), then another cut starting from the threshold of the 
three-pion production, that contributes to the isoscalar part of the nucleon EE's (in a VMD framework, 
one has the production of a virtual u meson, that decays in three pions), etc. (see below). In summary. 



within a DR approach one obtains (cf. Eqs. (39) and 
EE's, valid for any 



(41)) the following expressions for the nucleon 
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IT 



Ami. 



s{s — q^ ~ It) 
^m{F^^{s)] 
s{s — — ie) 



(48) 



where the normalizations of the EE's have been explicitly shown, with 



IS _ l^P~^ f^r. 



IV f^P f^^- 



Moreover, the lowest thresholds for each isospin channel have been considered. Three comments are in 
order: i) the main task should be the evaluation of the imaginary part of the form factors for s > 4m^, 
ii) one could expect more cuts starting at various thresholds, depending upon the particle productions 
one has to consider, iii) one should keep in mind that the dispersion integrals have contributions also 
from values of s below the NN production threshold (4m^ < s < 4M^), the so-called unphysical region. 

Indeed, the analysis of the discontinuity across the real axis, for 3f?e > 4m^, of the current-operator 
matrix elements allows one to determine the imaginary part of the EE's, establishing the announced 
link with physical processes, like the hadronic decays of neutral vector mesons, the TL form factors 
of mesons (e.g. for describing the 7r"'"7r~-continuum) or the physical NN annihilation into mesons. To 
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show this, still following Ref. one can write 



I EE' 



i^nf {0\Jif^iOmP,sWip',s')) 



E 



M. 



N 



{2nf' I / d'x e-'^'-^ vAv'n-t) (0| &(0), Sm{x)] W{p, s)) 
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vAp') {0\[JLiO),S^ix)] \Afip,s)) 



E 



(27r 



,3/2_ 



(2vr) 



du d'x e-"^P'+P-P^>'' 



-iujt 



u — te 



X vAp') {0\JUO)\pr,I,B = 0){pr,I,B = 0\SN{0)mP^s)) 

^3 ^ 



— (27r)3/2 
Mjv 27r 



P/-P -pj 



El - E' - E - ie 

X (0| J,';„(0)|p,; J, 5 = 0)(p,; /, 5 = 0|5^(0)|Ar(p, s)) = 

V 



Vs'ip) 



— (27r)3/2 
Miv 27r 



P/-P -PJ 



Er-E' -E 



+ i7i6{Ei -E' -E) 



X vAp') (0|J,^^(0)b,;/,i? = 0)(p,;/,i? = 0|5jv(0)|Ar(p, s)) (49) 

Notice that i) Eq. ( [46| has been used, and ii) 5'Ar(0) lowers the baryon number by a unity and therefore 
the intermediate states must have -B = 0, as also requested by the matrix elements of the em current 
(0| J^^(0)|p/; /, 5 = 0). If one repeats analogous steps (illustrated in Appendix [A|) for the matrix ele- 
ments of (0| [J^m(O)]^ W{P-, s)N'{p' ■, s')) = {M{p, s)J\f{p', s')| J^^(0)|0)*, one can single out the imaginary 
part of the nucleon FF's by evaluating matrix elements of the skew-hermitian part of </^(0), as follows 



I EE' 1 

^{2^f - |(0|JL(0)|Ar(p,.)Ar(p',.')) - (0| [JL(0)]^ \M{p,smp\s')) 



E 

Ml 

Vs'ip') 



Us{p) (50) 



From Eq. (50), by applying the proper Dirac traces, one can obtain the needed relations for determining 



53m |f^2)(5'^)| (see the next Subsection). The attractive feature of Eq. (50) is the link with physical 

processes, as already said. In particular, the matrix elements (0| J^(0jlp7; /, -B = 0) is related to 
the em production of physical hadronic systems with B = 0, while the matrix elements {pf, I,B = 
0\Sn{0)\N'{p, s)) is proportional to the hadronic production of a NN pair (or the NN annihilation into 
hadronic states). It should be pointed out that, since the intermediate states have pj = {p + p'Y > 0, 
given the presence of the four-momentum conserving delta-function, they represent physical states 
(differently from the SL case where {p' — pY < 0, i.e. unphysical intermediate states), and a spectral 
decomposition of the nucleon FF's can be devised. Strictly related to the value of the masses of the 
intermediate states, there is the question of the energy thresholds of the multi-cut structure that must 
appear in Eq. (48). One can immediately realize that the em hadron production, (0| J^„(0)|p/; /, 5 = 



0), has 2m7r as the lowest physical threshold (contributing to the isovector FF's), while the physical 
threshold of the hadronic production of the NN pair (p/; I,B = 0\S]\f{0)\M{p, s)), is 2Miy. The region 
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between these two thresholds, as mentioned above, is the unphysical region, that, however, contributes 
to the dispersion integral, since the path of integration starts at 2m^, i.e. where the free propagation 
of the pions begins. Therefore, one of the issues in the DR method is given by the evaluation of the 
matrix elements {pi;I,B = 0|S'Ar(0)|A/'(p, s)) in the unphysical region. This can be overcome through 
an analytic continuation of the physical processes and by considering the crossing symmetry (see |1] 
and |108] for a more recent discussion). Below the threshold 2Mj\f, there are other possible energy 
thresholds related to the em productions of three pions, i^^-pair, four pions, KKn, etc. In particular, 
one has to remind that the intermediate states, \pi',I), must have vanishing baryonic number and 
strangeness, and therefore they can contain two or more mesons, with proper quantum numbers, or 
baryon-antibaryon pairs. Moreover, given the presence of the matrix elements of the current operator 
(0| J^^(0)|p/; I , B = 0, S = 0) , the intermediate states must have the proper isospin quantum numbers 
dictated by Jemi^)^ that contains an isoscalar term and an isovector one. For instance, the first term 
filters non-strange meson states with I^[J^^) = 0~(1 ), i.e states with an odd number of pions, while 
the second term can single out non-strange meson states with J'^(J^'^) = 1^(1 ), i.e states with an 
even number of pions. Since the G-parity does not affect states made by strange mesons, the lightest 
state with a kaon and an antikaon can contribute to both isospin channels. 



3.6 TL Nucleon form factors and the dispersion relation approach 

The DR approach has been widely adopted for describing the SL nucleon FF's (see, e.g., |1U9] and 
references quoted therein). Particular care was devoted to modeling the spectral representation of the 



imaginary part of the FF's (cf. the sum in Eq. (50)). One of the main ingredients for achieving such 



a description was represented by the phenomenological input supplied by the VMD (see, e.g., [3]), 
since one could use a physically meaningful scheme and exploit experimental quantities for describing 
'^m{F^2){(l'^)} ■ Moreover, the extension to the TL FF's was straightforward, since the DR's hold for 
any value of g^. In particular, the DR approach adopted in Ref. |109j for a very accurate description 
of the SL FF's (as known at the date), was extended in Ref. |11U] to the TL region, without changing 
the number of the free parameters. For brevity we call this approach the Mainz model. In the spectral 
decomposition of 9'm{F^2)(5'^)}, it was considered: i) the same numbers of VM poles as in |109] . 
ii) the 27r continuum, already used in the SL region for constraining the overall description through 



unitarity, and iii) the asymptotic behavior suggested by pQCD (cf. Eq. (12)), including in this way both 
hadron and quark degrees of freedom. For the sake of concreteness, let us mention that the spectral 
decomposition in terms of the VMD model amounts to a finite sum on the intermediate states with 
masses equal to the relevant vector mesons, viz 

^mF^MD-At') = ^ E <n S{t' - {mif) (51) 



where i = 1,2 and / = IS, IV, and af „ contains all the information about 7* — VM coupling and the 
transition VM — )■ NN. The above result follows from the VMD contribution to the whole FF's, viz 



I 



where t = q^. Notice that the denominators, that come from the VM propagators, vanish only for t > 0, 



and produce the well-known pole pattern, experimentally seen in the TL region(cf Subsect. 2.2). 
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In |llU j. it was assumed an unsubtracted form for DR's, that allows one to write 
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where the function given by 
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(53) 



(54) 



takes into account the pQCD prescription for the asymptotic behavior. In the second line of Eq. (53) 



Ff'{t) describes the two-pion contribution, including both the term from the p pole and the 27r contin- 
uum. This last term is parametrized through the pion form factor and the P-wave tttt — )■ NN partial 
wave amplitude, without free parameters (see |109j for details). It is worth recalling that by using 
Eq. (54) one can recover the standard pQCD behavior for t — t- — oo, viz 



^e{Fi{t)} ^ 



In 



\t\ 

Ql 



^ -7 



(55) 



that, in turn, can be properly continued in the TL region. In Eq (54), 7 is the anomalous dimension, A 
controls the transition to the pQCD behavior and Qo ~ ^qcd- In the actual calculations, the constants 
were chosen as follows: 7 = 2.148, corresponding to a number of flavors Nf = 3, = 0.35 GeV^, 
while the adjusted parameters were: i) the quantity A and ii) the residues af„, corresponding to 
the contributions from three meson poles for each isospin channel. In particular the masses of the 
vector mesons were: m(p') = 1.40 GeV, m(p") = 1.45 GeV, m{p"') = 1.69 GeV, m{uj) = 0.782 GeV, 
m{(j)) = 1.019 GeV and m(S") = 1.60 GeV. Results for the TL nucleon FF's, obtained by applying the 



Mainz model, are shown in Fig. (27) 



The work described in Ref. [llOj has been steadily improved by the Mainz-Bonn- Jiilich collabora- 
tion |llllll08j . enriching the previous description both i) by adding the KK as well as the pn continuum, 
and ii) by updating all the experimental inputs with latest experimental results. Moreover, different 
ways for implementing the asymptotic behavior requested by pQCD (see Refs. |llll I108j for more de- 
tails) were considered. In particular, in the most recent SL and TL calculations presented in Ref. |108] . 
it was introduced the so-called superconvergent approach for controlling the asymptotic behavior of the 



TL nucleon FF's, besides the standard one given by Eq. (55). The superconvergent approach amounts 
to impose a proper fall-off not only to the real part (cf. Eq. (55)), but also to the imaginary part of 
Dirac and Pauli FF's, i.e. 



POO 

/ dt Qm{Fi{t)} = 

Jto 



dt 3m{F2(t)} = 



to 



dt t 3m{F2(t)} = 



(56) 



to 



It should be pointed out that these constraints allow to take into account the application of the 
Phragmen-Lindeloff theorem, as illustrated in Subsect. 3^ The TL nucleon FF's, calculated by 
using such a superconvergent approach are shown in Fig. (28). Notice that i) the SL and TL proton 
data and ii) the neutron SL FF's have been included in a 17-parameter fit, and therefore the neutron 
TL calculation (solid curve in the rightmost panel of Fig. (28)) represents a genuine prediction. In 
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Figure 27: TL nucleoli form factors within the Mainz model |110j . Solid lines: calculations with three 
meson poles, for each isospin channel, and = 12.0 GeV^ (see text). Dashed lines: the same as the 
solid lines, but including the contribution from 0(1680) as well. Notice that the experimental data 
are presented as the absolute value of magnetic form factors, but the experimental data correspond to 
effective proton form factor given by Eq. Figure taken from Ref. |110j. 



Ref. |108] . it was also presented an explicit pQCD continuum approach, based on a phenomenological 
Ansatz that fulfills the pQCD behavior (55). In this case the set of free parameters is reduced to 14, 
but the quality of the fit become worst, as pointed out in the same work. 



More recently, the above illustrated approach has been improved by the Tokyo-Chiba group, adding 
more vector-meson poles and more parameters for modeling the pQCD contribution. In Ref. |112j . an 
unsubtracted DR has been assumed for both Dirac and Pauli form factors, in the SL and TL regions. 
The spectral decomposition of the '^m{F^^^^{q'^) is described in an analogous way as in Ref. [TDSJ. It 
was considered i) a 27r continuum with a helicity decomposition of the vrvr — > A^A^ amplitude (notice 
that the tabulated data for the vrA^ scattering adopted in |112] . i.e. the ones by Hoheler et al |113] . 
already contain the p contribution), and ii) up to 5 meson resonances for each isospin channel, taking into 
account both mass (up to mvM ~ 2.05 GeV) and width in the Breit-Wigner function that describes each 
vector-meson pole contribution. A particular care was devoted for modeling the asymptotic behavior 
suggested by pQCD, both considering contributions up to ^^(t) and implementing the superconvergent 
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Figure 28: TL nucleoli form factors evaluated within the Mainz-Bonn- Jiilich Super Convergent disper- 
sion relation approach |108] . Left panel: proton form factor in the range 3.5 < t < 5.5 GeV^. Middle 
panel: the same as the left panel, but in the range 5.5 < t < 20 GeV^. Right panel: neutron form factor. 
Solid lines: best fit. Dashed lines: error bands assigned by allowing some variation to the y^/DOF 
(see [lOSj for details). Notice that the experimental data are presented as the absolute value of magnetic 
form factors, but they corresponds to the effective nucleon FF, introduced in Eq. ([2]). Figure reprinted 
with permission from Ref. |108j . ©2007 American Physical Society. 



Ansatz in a slightly different way with respect to Ref. |108] (cf. Eq. (56)), viz 

- / dt' {t'Y "^m {Flit')} =0 n = 0, 1 
vr 74^2 

- / dt' [t'Y 53m{F/(t')} =0 
where Q^m |F/(t)} has to satisfy the asymptotic condition 



n = 0, 1,2 



(57) 



^m{Fl{t)]^ 



In 
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with c and 7 constants. In the actual calculation it was used 7 = 2. With the previous constraints, 
Fi{t) has the standard pQCD behavior, a part the value of 7 different from the standard anomalous 
dimension. In the approach of Ref. |112] , the number of free parameters has substantially increased, 
since one has 5 masses and widths for each isospin channel, and correspondingly 10 residues for the 
vector meson poles of both the Dirac and Pauli FF's. Therefore the VMD part contains 40 parameters. 
Finally, one has to add 12 adjusted coefficients for the QCD part, 3 parameters that enter a convergence 
factor, ensuring the super-convergence condition and another one that shifts slightly the masses of the 
resonances in the Breit-Wigner width. The 56 parameters are fixed by fitting: i) the most recent SL 
and TL proton data, even considering the ratio jG^/G^jl measured by the BaBar collaboration [Qj in 
the TL region (notice a mistyping in |112] . where the experimental ratio is indicated with /XplG^/G^I) 
and ii) the SL and TL neutron data. The reasonably good description of the nucleon FF's, achieved by 



the fitting procedure of Ref. |112j . is illustrated in Figs. (29) and (30). 
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To summarize the investigations presented in this Subsection, one has to stress the richness of 
information in the TL region, that can be accessed through Qm{F^2)} . Noteworthy, the experimental 
determination of the relevant phases could strongly motivate future efforts in developing more and more 
refined theoretical analysis, possibly tightly related to microscopic calculations, in order to establish a 
close vision on the non-perturbative regime of QCD. 
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Figure 29: TL proton effective form factor, Eq. (25), within the fitting procedure of Ref. |112] . based on 



the DR approach and a super-convergence condition (cf. Eq. (57)), vs the Mandelstam variable s 
Left panel: < 5 GeV^. Right panel: the whole range of s. Figure reprinted with permission from 
Ref. [112], ©2010 American Physical Society. 
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Figure 30: The same as in Fig. (|29j), but for the neutron effective form factor, Eq. (25). 
reprinted with permission from Ref. |112j . (c)2010 American Physical Society. 



Figure 



3.7 Phenomenological DR's and TL extrema 



In Ref. |114j the nucleon magnetic FF's below the A^A^ threshold, i.e. in the region < g 
has been investigated in an almost model-independent way. The DR's are adopted for determining 
an integral equation that relates log Gn{(i'^) and log\GN{q'^)\ where GN^q^) = Gui^iq^) I IJ'N (therefore 
GAr(O) = 1), following the modulus representation of the FF discussed in Ref. |115j . In the TL region 
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one writes the following subtracted (recalling that log Gn{0) = 0) dispersion relation 



Zo^ Gjv(g ) = log \GN{q )\ + ?$G^(g ) = — / dt — , (59) 

where $G^(g^) is the phase of GM{q^). The range of integration can be split into two parts, 4m^ < 
t < AMjq and 4M^ < t < oo. The integral in the second interval, i.e. the physical region, has been 
estimated by adopting a suitable function for \GN{q'^)\, with the expected asymptotic behavior, and 



fitting the available experimental data, assumed to directly give \G^{q'^)\ (cf. Eq. (25) for the proper 
caveat). Then the value of log G]y{q'^) in the unphysical interval, < < 4M^, is obtained by solving 
the following integral equation, through a very detailed procedure in order to overcome many numerical 
subtleties, 

<o« G.if) - IM) - '^^^^ r * . , (60) 



TT J 4ml t (t-g2)^t-4m2 

where Iph{,q^) indicates the above mentioned integral that depends upon \GM{q^) \ in the physical region. 
Interesting structures, related to peaks for i) p(770), ii) p'(1600) and iii) near the NN threshold are 
automatically generated. Noteworthy, no evidence of a bump at the 0(1020) mass has been found, in 
spite of the expectation that such a peak should exist if a sizable polarized ss content in the nucleon 
were present. 

Within the DR spirit, in Refs. |116t 1117] . the proton ratio Rp{q^) = IJ^pG Ep{,q^) / G Mj,{,q^) has been 
investigated by considering the following dispersion relation 

R^if) - n^iO) + ^ r * (61) 

TT J4^2 t [t - q^) 

where Rp{0) = 1. Moreover, the Phragmen-Lindeloff theorem was implemented by imposing a vanishing 
asymptotic value for Q'm{i?p(g^)}, in order to recover a real function when — )• oo (as in the SL region). 
As a final step, the imaginary part has been spanned onto the ortho-normal basis given by the Chebyshev 
polynomials, with independent coefficients for the unphysical and the physical regions. The two set of 
coefficients are fixed by means of the available experimental data and the theoretical constraints in 
both the SL and the TL |116[ I117j . obtaining a description in the unphysical region and also in a wide 
range of g^. Unfortunately, the fitting procedure cannot be fully constrained, and this is refiected in the 



wide band associated to the solution of the integral equation (61). However, interesting results have 
been obtained, both near the threshold and in the asymptotic region. In particular, a dominance of 
Gg(g^) near threshold has been found, as indicated by a growing ratio up to ~ 11 for g^ ~ 4 (GeV/c)^. 
As to the asymptotic region, the moduli of the ratio are quite similar, but the prediction based on 
the Phragmen-Lindeloff theorem cannot be immediately tested. As a matter of fact, the SL result for 
Rp{Q^) is negative for Q"^ = — g^ — )■ oo, and in the TL region the extracted value is only the modulus 
\Rp{q'^)\. Therefore the phase ^(cx)) of Rp{q^) = |-Rp(g^)| exp[i$(oo)] for g^ — )■ oo has to be determined, 
for a full comparison. Indeed, in order to fulfill the equality between the SL and TL asymptotic 
values, i.e. limq2_>.oo 3fJe{-Rp(g^)} = 3fte{-Rp(— g^)} (as dictated by the application of the Phragmen- 



Lindeloff theorem; cf. Subsect. 3.4), one should expect that $(oo) = {2n + 1) vr, since 9fJe{-Rp(— g^)} 
asymptotically is negative. From such an observation, one can deduce the following interesting result. If 
one applies the argument theorem |105j to a closed path, composed by a circle at infinity and the upper 
and lower rims of the cut [4m^,oo], and reminds also the Schwarz refiection principle, one gets that 
$(oo) — $(4m^) = 7r(A^ — P), where N{P) is the number of zeros (poles) inside the closed path. Notice 
that P = 0, given the assumed analyticity inside the closed path and one can choose $(4m^) = 0. 
Therefore, N must be an odd number to fit a negative value for 3fJe{-Rp(g^)} = |-Rp(g^)| cos[<l>(oo)]. 
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This seems to be in nice agreement with the SL data suggesting the existence of at least one zero (cf. 



Rp{Q ) in Fig. (23)). Further analyses of the issue have been carried out in Ref. [11 7j . where it has 



been discussed the large uncertainty that affects the modulus of the previous ratio, due to the choice of 
the data set that one adopts for constraining the high momentum part of '^m{Rp{q'^)} . In particular, 
one gets 



0.95 ±0.20 
2.3 ±0.7 



from BaBar [9] 
from Lear ITOl 



Summarizing, all the previous results indicate that the phenomenological studies performed through 
the DR's are a source of highly non trivial information in both the extrema of the TL region, given the 
powerful tool represented by the complex analysis. 



3.8 The analytic continuation of the SL Nucleon FF's 

In the DR approach, in its full glory, the main issue is given by the attempt of exploiting physical 
processes for modeling the imaginary part of the nucleon FF's, or better the discontinuity of the cur- 
rent matrix elements across the multi-cut region for 3fJe > 4m^. As we have seen in the previous 
Subsections, the knowledge of the imaginary part of the nucleon FF's allows one to determine the 
FF's themselves in the whole range of g^. A different program is presented in Refs. |118[ 1119] , where 
it was proposed a fit of the TL data, starting from a well-known model of the SL FF's, like the one 
in Ref. [66], devised in the 1970's for describing, through a VMD approach (see, e.g., Ref. [3]), the 
available SL nucleon FF's, but able to predict the fall-off of the SL ratio UpG^^iQ"^) /G^jiQ"^) recently 



observed (cf. Subsect. 3.1 and 0). The simple strategy for generating the TL nucleon FF's, that are 
complex functions, is based on the analytical continuation of the SL FF's, that are real functions of 
3fte q^ < 0, replacing the dependence upon = —q^, with —q^ exp{i9) and then using the phase 6* as a 
free parameter in the fitting procedure to obtain the TL nucleon FF's. VMD plays a fundamental role 
in the DR approach for reconstructing $5mFi(2)(g^), since it provides the pole structure, to which one 
has to add the multi-cuts generated by the free propagation of the hadronic states. Therefore the VMD 
content of the SL FF's allows one a first approximation, that can be enriched by using the adjusted 
parameter 6. Indeed, only an infinitely accurate knowledge of the SL FF's (i.e. only in a mathematical 
realm) could univocally determine the FF's in their analytical domain, but this is not the case within 
the realm of the phenomenology. 

Refs. I119j have generalized the SL model of Ref. [66j, where the external photon couples to 
both an intrinsic structure, described by a proper form factor and a meson cloud, obtained within the 
VMD framework. It should be pointed out that the presence of the intrinsic form factor has a clear 
physical motivation: the size of the constituent quarks (CQ's) inside the nucleon (see, e.g., Ref. [120[ 
for a discussion of the CQ FF in a relativistic description of the SL FF's with a realistic nucleon 
wave- function) . In particular, the SL model is given by 
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Figure 31: TL nucleoli FF's calculated within the approach of lachello et al. |118[ I119j (see text). 
Upper panel: the magnetic form factors \GMp{(f')\ and \GM„{'f)\- Lower panel: the electric form 
factors \Gep{(1^)\ and \GE„{(f)\- The solid lines are from the analysis of Ref. [119] and the dashed 
ones are from previous calculations presented in Ref. |118] . The experimental data correspond to the 
effective TL nucleon FF's, as given in Eq. ([2]). Figure reprinted with permission from Ref. |119] . ©2004 
American Physical Society. 



where 



(1-7 q^Y 



is the intrinsic CQ FF, nip = 0.776 GeV, rrica = 0.783 GeV, = 1.019 GeV. The six free parameters, 
/3p, ap, and 7 are fitted to the SL data. It is worth noting that the asymptotic behavior of the 

previous parametrization is the same predicted by pQCD, without logarithmic corrections. Moreover, 
in order to take into account the non negligible width of the p meson, the propagator is changed 
accordingly to 

rUp m^ + 8rpm^/7r 



m2 - g2 ^2 _ q,2 _j_ ^4^2 _ q'i'^Y p[a{q'^) - il3{q'^)]/m^ 

with a(g^) and being adhoc functions. The value of the fitted phase was found to be 6* ~ 53° in 

Ref. |118] . and ~ 22.7° in Ref. |119] . where the proton data at the largest values were not included. 
In the upper panel of Fig. (31), the theoretical magnetic FF's, \GMp{,q^)\ and |G'm„(q'^)|, are presented 
and compared with the experimental data, that correspond to the effective nucleon FF's (see Eq. (25) 
and the relative caveat). In the lower panel, the TL electric FF's are also shown. Unfortunately, as 
illustrated in the case of \GMp{q^)\, the model does not contain the structure in resonances shown by 
the TL data for q^ > 4M^ (the adopted VMD stops at the 0- meson), and therefore the comparison 
with the experimental data is not satisfactory. As to the neutron, it should be noted the large value 
reached at the threshold, in agreement with the FENICE measurements [T3] . 

In the same spirit as the previous approach, in Ref. [70] the TL generalizations of both the SL model 
of Ref. [66] and the SL model of Ref. [67] have been elaborated. In the last case, the 14 parameters 
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Figure 32: TL nucleoli FF's evaluated within the approach of Tomasi-Gustafsson et al. ^Q\. Dashed 



lines: pQCD-inspired model given in Eq. (64). Dotted lines: TL extension of the SL VMD model of 

j . Notice that the experimental 



Ref. [66] • Solid lines: as the dotted lines, but for the SL model of Ref. 
data corresponding to the effective FF given in Eq. (|2]), are separately compared to the magnetic and 
electric Sachs FF's. Figure reprinted from Ref. [7D] with kind permission of The European Physical 
Journal (EPJ). 



present in the SL nucleon FF's have been updated by fitting the most recent SL and TL nucleon FF's 
data-base (see Ref. |i70j for the full list). In particular, the generalization to the TL region was performed 
exploiting the following relations 



Q==-«==«'e--»(";|j^"i;i:;j (63) 



Together with the above mentioned models, the analytical continuation of a pQCD-inspired model was 
also proposed, based on counting rules and following the dipole behavior in the SL region. Such a model 
becomes in the TL region 

where A = 0.3 GeV is the QCD scale parameter, A{p) = 56.3 GeV^ and A{n) = 77.15 GeV^ are 
extracted from the fitting procedure to the TL data. The results of the approach of Ref. [70] are pre- 



sented in Fig. (32 ). Notice that the SL model of Ref. [67] contains two extra vector mesons with masses: 
nip' = 1.45 GeV and m^^/ = 1.419 GeV. Unfortunately, also in this case, the TL data, corresponding to 
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the effective nucleon FF's, given by Eq. (|2j), have been compared with both \Glj\ and \G^\, separately. 
In view of this, the comparison proposed in Fig. (32) should be carefully considered at large values of 

Very recently, within the analytic continuation framework, an updated version of the FF's of Ref. [6 
have been applied in Ref. |121j for obtaining a unified fit of both SL and TL FF's. 



3.9 Nucleon FF's within a Relativistic Constituent Quark Model 

A more microscopic description of the nucleon FF's, in both SL and TL regions, was proposed in Ref. [7], 
where it was adopted a relativistic Constituent Quark Model (RCQM), already successfully applied to 
the pion FF in the whole kinematical range |B]|. The aim of the approach is to calculate as many 
quantities as possible in terms of quark degrees of freedom and to perform a direct evaluation of the 
SL and TL FF's, without taking care of the DR constraints. At the end, after fixing only four adjusted 
parameters by using the SL experimental data, one can predict the proton and neutron FF's in the TL 
region. In Ref. [7] , the following ingredients were chosen: i) the Mandelstam formula |122j for calculating 
the matrix elements of the current operator in terms of the nucleon Bethe-Salpeter amplitude, ii) a 
microscopic relativistic Hamiltonian that is able to describe reasonably well the isovector and isoscalar 
vector meson spectra |123] . to be included in the VMD part of the quark-photon vertex, iii) a quark- 
photon vertex decomposed in a bare term and a VMD term, iv) a Light-front (LF) framework |125] . 
for embedding the previous ingredients and obtaining a relativistic (LF-boost invariant) CQM. 

In the RCQM of Ref. , the quark- nucleon vertex is described through the Bethe-Salpeter ampli- 
tude, with a properly symmetrized Dirac structure, approximated by 

^■'{h,k2,h,p) = I [ S{h) Ty 7^ Scih)C ® Sih) + S{h) Ty 7^ Sc{h)C ® Sih) + 

+ S{h) Ty 7^ Sc{k2)C ® Sih)] Aih, k2, h) xr^ usip) (65) 

where Xr^ is the isospin spinor, S{ki) the Dirac propagator of the quark, Sc{ki) = C S'^iki) C^^ = 
75S'(/cj)75 with C the charge conjugation and A(A;i, ^2, ^s) yields the symmetrized dependence of the 
vertex function upon the quark momenta, ki. 

In the TL region (for the SL expression see Ref. [7]), the matrix elements of the macroscopic current 
can be evaluated microscopically by means of the Mandelstam formula jl22] . i.e. 

X 1 ^ 1 ^Tr{5^i(A:2)$?;(fc„fc2,A:^,p')^-^(fci)X'^(fc3,g)<f^(fci,A;2,fc3,-p)} (66) 

where in $'^(/ci, k2, k^, —p) appears Us{—p) = 'i^s(p), l-^ih-, q) is the quark-photon vertex, the number 
of colors and the trace runs over isospin and Dirac indexes. Notice that ki + k2 + k3 = p and ki + k2 + k'^ = 
-p'. 

The quark-photon vertex is decomposed as follows 

— j-jg -r I z-i-fv ' 

where is the third component of the isospin operator and each term has a purely valence contribution, 
acting in the SL region only, and a contribution corresponding to the em gg-pair production (non valence 
contribution). In turn, the pair-production contribution, i.e. the relevant one in the TL region, can be 
decomposed in a bare term plus a VMD term (according to the photon-state expansion given in Eq. 



(19)). For instance, in the SL region, for a CQ with momentum k'^ one has 

I^{k, q) = Afj 9{p+ - k^)e{k^)^^ + ^(g+ + k^)e{~k^) [Zb Afn^ + 4m ^tik, q) } 



valence non valence (67) 
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where = + k^ is the so-called plus component of the quark four-momentum (see Ref. [125J for a 
detailed review of the LF framework), / = IS, IV, Afis = 1/6 and Afjv = 1/2. It has to be pointed out 
that the constraints 6{k~^) and 6{—k^) single out a quark or an antiquark, respectively. The constants 
in the non valence term, Zb (bare term) and Zyj^j (VMD term), are unknown weights determined from 
a fit to the SL data-base (indeed, in the actual fitting procedure Zb = Zy)^^ was assumed, see below). 



In Eq. (67), T^j{k,q) yields the VMD contribution to the quark-photon vertex and it has the same 
structure of the one adopted in the pion case [6J, but including the isoscalar vector mesons. For a given 



isospin channel, T^j{k,q) is approximated as pictorially illustrated in Fig. (33) (cf. Ref. |i6j for more 
details in the pion case), viz 



[6^(P/;n)]7/;n 



g2 -m^.^ + z m7.„r7.„(g2) 



(6J 



where e\{Pi.n) is the VM polarization, is the em decay constant of the ra-th vector meson, m/.„ 
the corresponding mass, A/.„(fc,P/.„) gives the momentum dependence of the VM vertex function and 
Vj.^ik, k — Pj-n) is the corresponding Dirac structure, while ^/.^(g^) is the running total decay width. It 



has to point out that in Eq. (68) all the quantities appearing in the numerator have been approximated 
through the corresponding values calculated at the respective pole (i.e. the residue at 



m 



I:n 



obtaining, therefore, the correct result when the contribution is maximum. In turn, such an approxima- 



tion allows one to evaluate 



ex{Pn)-Vn{k,k-Pn) 



An{k, Pn), (our Ansatz for the VM vertex function. 

I.e. the amputated VM Bethe-Salpeter amplitude), by using the eigenvalues and the eigenvectors of the 
mass operator for vector mesons, proposed in Ref. |123] . 

The running total decay width in the denominator of Eq. (68), r7.„(g^), is vanishing in the SL region, 
while in the TL region it is assumed to be 



f/;n(g' 



I:n 



q^ — Ami 


3/2 




1/2 


.„ - Ami 




q2 





(69) 



where F/.^ is the total decay width (see, e.g., PDG |124] ). 

Once the on-mass-shell approximation has been introduced in Eq. (68), F^(A;,g) is not anymore a 
four vector, but since only the plus component Tf{k,q) = T^{k,q) + Tj{k,q) has to be used in the 
evaluation of the nucleon FF's, it is possible to generate matrix elements of the current operator that 
correctly transform under the LF-boosts (these boosts act on the plus component of a given four-vector 
without mixing other components, differently from the standard boosts) . 

Up to 20 isoscalar and isovector VM's have been considered, since in the case of the pion that 
number of isoscalar mesons was necessary in order to get fully convergent results for the TL pion FF 
at the largest experimental values of q^ {q^ ~ 10 (GeV/c)^). The first four (three) meson masses 
for the isovector (isoscalar) channel are taken from PDG |124j together with the corresponding total 
decay widths F/.^, while all the other masses come from the solutions of the mass eigenequation |123] . 
and F/.„ = 0.150 GeV is taken, as in the pion case. Noteworthy, all the em decay constants fi-n are 
calculated from the corresponding eigenfunctions, and they allow the evaluation of the em decay widths 
V^^+^_, that result in a very nice agreement with the experimental values, when available. Finally, as in 
the case of the pion, the quark mass mq = 200 MeV is adopted. In conclusion, the VMD part of the 
photon-quark vertex is fully constrained, a part the relative weights Zy^j^j and Zy^^j = Zb- 



The momentum dependence of the nucleon vertex function in Eq. (65) is differently modeled in 
the valence region (necessary for evaluating the SL nucleon FF's) and in the non valence one. In the 
valence sector, the spectator quarks are on their-own mass-shell, and the 3-momentum dependence, is 
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Figure 33: VMD contribution to the quark-photon vertex. The double-wiggly hues represent the Green 
function describing the propagation of the n-th VM. The loop on the right indicates the decay constant 
of the n-th VM. The LF-time flows from the right toward the left. Figure reprinted with permission 
from Ref. |6], ©2006 American Physical Society. 



approximated through a nucleon wave function a la Brodsky-Lepage (pQCD inspired) |125] . given by 



[M^-M2(l,2,3)] 



{iiUzY [/32 + M2(l,2,3)] 



7/2 



(70) 



where Mo(l,2,3) is the free mass of the three-quark system, = kf / and M a normalization 
constant, fixed through the normalization of the proton charge FF. The power 7/2 and the parameter 
p = 0.13 are chosen for obtaining an asymptotic decrease of the valence contribution to the SL nucleon 
FF's faster than the dipole. The value (3 = 0.65 GeV can be fixed through the evaluation of the 
anomalous magnetic moments, fip and fin, obtaining for those quantities, fip = 2.87 (Exp.= 2.793) and 
/i„ = -1.85 (Exp.= -1.913). 

The non valence vertex, that is fundamental for obtaining a zero in the SL ratio ^pG^^{Q'^) / G\,j{Q^) 
(cf. Fig. (23)) and essential for evaluating the TL FF's, can depend on the available invariants, e.g. in 



the SL region one has i) the free mass of the spectator quarks, Mo(l,2), and ii) the free mass of the 
N-g system Mq{N, 3). Then, in the SL region, the non valence vertex is approximated by 



a5L 



{ki, k2, ks) 



[guf [gm 



|7/2-2 



k 



12 





r 











(71) 



where /c^ = k^ + /cjand qab = (i^a ^b)/ + ^o(^' -^)]- "^^^ power r is a free parameter. In the 
TL region, the non valence vertex can depend on the mass of the nucleon-gg system, Mq{N, 12), and it 
is given by 



xTL 



^[gi,2ng.T2Y^'-' 



-k 



12 



N 



P 



N 



k' 



P. 



N 



k' 



(72) 



where the factor of 2 counts the possible patterns for the gluon emission by the produced qq pair. 

The free parameters determined through a fit to the available SL experimental data, G\,j{Q'^), 
G^i^Q"^), fipG^^{Q^) / GjAQ'^) and G^(Q^), are four: i) two weights out three present in the pair pro- 
duction term of Eq. (67), i.e. Zb = Zy^ = 2.283 and Z(f^/Z(y^ = 1.12; ii) the power p = 0.13 in 
the valence amplitude, Eq. (70); iii) the power r = 0.17 in the SL and TL non valence vertexes, see 
e.g. Eqs. (71) and (72). The minimization yields = The overall quality of the SL fit can be 



appreciated from the comparisons in Fig. (23), and from both the calculated proton charge radius Vp 



(0.903±0.004) fm {rf'^ = (0.895±0.018) fm), and the neutron charge slope at q^ = 0, -dG%{q^)/dq^ = 
(0.501 ±0.002) (GeV/c)-2 (exp value; -dG%{q^) / dq^ = (0.512 ± 0.013) (GeV/c)-^). Surprisingly, the 
same quality in the fitting procedure can be achieved by excluding the UpG^^iQ"^) / G^^^iQ"^) data [7]. 
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Figure 34: Experimental nucleon effective form factors, Eq. ([2]), in the TL region compared to the 
predictions of the relativistic constituent quark model of Ref. [7], evaluated according to Eq. (25). 



Left panel: proton FF (data from [9] and references quoted therein); notice the logarithmic scale for 
q^. Right panel: neutron FF (data from Ref. jl3]). Solid hne: RCQM full calculation. Dashed hne: 
neutron full calculation multiplied by an arbitrary factor of two. Figure taken from Ref. [7]. 



In Fig. (34), the experimental effective FF's (see Eq. rt2])) for both proton and neutron are compared 
with the RCQM results (evaluated according to Eq. (25)), that are true predictions of the model. It is 
worth noting that a microscopic calculation could lead to the extraction of valuable information on the 
dynamical content of a realistic RCQM. For instance, missing strengths at = 4.5 (GeV/c)^ (nicely, 
the same feature appears in the TL pion FF [U]) and = 8 (GeV/c)^, could point to an enrichment of 
the present VM data-base, and therefore to a strong request of refining the RCQM approach in the VM 
sector. This simple example sheds light on the potentially relevant results that one could achieve from 
the forthcoming accurate measurements in the TL region and more and more improved microscopic 
calculations. 
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4 Conclusions sand Perspectives 



This review has illustrated the advancements achieved so far, both experimentally and theoretically, in 
the investigation of the timelike electromagnetic form factors of the nucleon. We have also presented 
the forthcoming experimental facilities, from which in the near future results are to be expected. 

The development of the field, started with the pioneering experimental work carried out in Fras- 
cati [m [13] up to the forthcoming exciting experiments with high intensity e^e~ beams ES] as 
well as - hopefully polarized - antiproton beams |62l ES] at FAIR, is demonstrating the continuous 
interest in these challenging studies. However, to measure separately the electric and magnetic nucleon 
FF's in the TL region is not an easy task, and the present-day database is still given in terms of an 



effective TL nucleon FF, that for the electron-positron annihilation process reads (see Eq. (25) and the 
corresponding remarks) : 



|Gf,(g^)| 



1 + ^ 



where the rightmost-hand side shows the relations with the Sachs FF's in one-photon approximation. 

The first success in the experimental investigations was the study of the total cross section for the 
reaction e^e~ — )■ by exploiting the energy scan technique at electron-positron colliders. In this 
respect, a major improvement was achieved by introducing the initial state radiation technique, i.e. by 



using events like e"'"e~ — )■ NN^y (cf. Subsect. 2.1.3). It was possible via ISR to obtain a measurement 
of the proton form factor with so far unprecedented precision and covering a very wide range. It 
should be pointed out that the ISR technique is of great appeal also for the planned high-luminosity 
e~^e~ B-factories, which are mainly devoted to the search of New-Physics signatures in B-decays. The 
study of the process cr'^+^^_^j^j^, lead to single out three energy regions: i) the region around the 
production threshold, ~ 4M^, that yields valuable information on the transition from unphysical to 



physical regions (cf. Subsects. 3.3 and 3.7); ii) the region where the Vector Meson Dominance and the 
interplay with the multi-cut analytic structure of the FF's fully display their phenomenological strength 
in the studies of the non-perturbative regime of QCD; iii) the asymptotic region where the realm of 
the perturbative QCD settles and the application of the Phragmen-Lindeloff theorem is representing 
a fundamental constraint for calculations based on microscopic models. As shown in the theoretical 
interpretation of the experimental effective FF, (let us recall that such an interpretation is driven by 
the one-photon approximation and that the two-photon corrections seem to be small in TL region 
[Ml UHl the electric and magnetic Sachs FF's are both present in a proper combination (cf. the 
rightmost side of the above expression). As previously mentioned, experimental efforts already have 
been carried out [101 E] and will be pursued in near future with the goal to fully disentangle the electric 
and the magnetic FF's. To this end, extending the measurements of the angular distributions (see [9| 
for recent experimental achievements for the proton case) could provide a tool for determining the ratio 
\G^/G^\. The holy grail in the field of TL nucleon FF's, however, is the investigation of the phases of the 



complex- valued FF's. This is heuristically illustrated in Fig. (35), where the discrimination power of an 
accurate measurement of the proton polarization perpendicular to the scattering plane, is shown. Notice 
that Py does not depend upon the polarization of the incoming electron beam, and once a polarimeter 
for the measurement of the proton polarization will be included into the experimental apparatus at an 
electron-positron collider, this measurement of Py could be performed. The experimental difficulties 
for achieving such a goal are known and are related to the measurement of the polarization degrees of 
freedom, but the strong motivations to pursue such a goal were already contained in the seminal work of 
Ref. [1], where the bridge between the TL FF's , in particular Qm{F^2)}> and the hadronic amplitudes 
of hadronic physical processes was established. Therefore, when measurements of 53m{F^2)} will be 
available, the possibility to extract accurate information on both the em decays of vector mesons and 
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Figure 35: Proton polarization, Vy (see Eq. (30)), normal to the scattering plane calculated by 
using different models. Dot-dashed line: the analytic continuation of the SL proton FF of lachello et 
al. [M| I118j . The other curves represent the analytic continuation of fits to the SL ratio Ff/F^, as 



indicated in the legend. 
Society. 



Figure reprinted with permission from Ref. |£8j, ©2004 American Physical 



the hadron production (with given quantum numbers) via NN should allow a carefully study of the 
hadronic states. From the theoretical side, the starting point of this program has been the analytic 
behavior of the complex-valued nucleon FF's, as imposed by the causality, and the consequent dispersion 
relations (cf. Subsect. 3.4), that link the real and imaginary parts of the nucleon FF's. The Lehmann- 
Symanzik-Zimmermann reduction formula made it possible, to express 53m{F^2)} in terms of the above 
mentioned hadronic amplitudes, that one can study in other experiments, since they are involved in 
physical (not virtual) processes. The first phenomenological attempts to model 53m{F-^2)}, as shown 
in Subsects. 3^, 3^ and 3^, were based on exploiting experimental inputs, like masses and widths of 
the vector mesons, but also i) introducing adjusted parameters for implementing a fine tuning and for 
describing the hadronic free propagation (i.e. the continuum associated to the multi-cut structure of the 
TL FF's) and ii) assuming theoretical constraints, particularly for the asymptotic behavior dictated by 
pQCD. The next step should be a refinement of our phenomenological tools, extending and enriching 
the microscopic calculations, like e.g. the one of Refs. [HI Ej, which are based on a LF Hamiltonian 
for describing the VM spectra and successfully applied to both the pion and the nucleon. Another 
appealing approach is the one based on the Anti de Sitter/Conformal Field Theory framework, where 
a very nice description of both SL and TL pion FF has been achieved [5] . In this way, by developing 
models with less and less phenomenological inputs, we could link the measurements of Q'mlF^g)} to 



55 



QCD-based calculations. 

In summary, a detailed comparison between accurate experimental data, in particular when both 
absolute values and phases of the TL nucleon FF's will become available, and microscopic approaches 
will allow to gain deep insights into the non-perturbative regime of QCD. 
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A Matrix element of the conjugated current 



In this Appendix, the steps for obtaining the equivalent expression given in Eq. (49), but for the matrix 
elements of (0| [J^rnWf \Af{p,s)M{p',s')) = {M{p, s)^{p' , s')\J^^{Q)\0)* are shown. One gets 



I EE' 



{2^f mp.s)u{p\s')\jumr 
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(27r) 



3/2 
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